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MUTUAL COUPLING IN FINITE 


ARRAYS OF RECTANGULAR APERTURES 

Abstract 

by 

JOHN WILLIAM SILVESTRO 

The mutual coupling between rectangular apertures in a finite 
antenna array, in an infinite ground plane, is analyzed using the 
vector potential approach. The method of moments is used to solve the 
equations that result from setting the tangential magnetic fields 
across each aperture equal. The approximation uses a set of vector 
potential modal functions to solve for equivalent magnetic currents. 

A computer program was written to carry out this analysis and the re- 
sulting currents were used to determine the co- and cross-polarized 
far zone radiation patterns. 

Numerical results for various arrays using several modes in the 
approximation are presented. Results for one and two aperture arrays 
are compared against published data to check on the agreement of this 
model with previous work. Computer derived results are also compared 
against experimental results to test the model's accuracy. These tests 
of the program's accuracy showed that it yields valid data. 

A look at the effects of mutual coupling in modest sized arrays 
and a short discussion on the beam scanning properties of three aper- 
ture arrays are also presented. Included with this data are results 
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for the first order coupling theory for small arrays. This theory 
assumes that the coupling of energy back to the source due to the re- 
radiation of nearby undriven apertures is small enough to be dropped 
from the analysis. The tests showed good agreements with the theory 
and are presented here for future use in the study of possible approx- 
imation techniques. 
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CHAPTER 1 


INTRODUCTION AND REVIEW OF PAST WORK 

1.1 Introduction 

It is projected that the growth of the satellite communication 
service will result in saturation of the existing frequency bands (C 
through Ku ) t by the end of this decade. In hopes of relieving this 
situation NASA is currently supporting research into the development 
of a 30/20 GHz satellite communication system. The aim of this 
work is to develop the technology necessary to make the use of these 
higher frequencies possible and at the same time to develop a more 
efficient system. In order to accomplish the latter, the satellite's 
antenna will have to provide narrow spot beams that are either fixed 
or scanned. With these spot beams, the antenna can service small 
areas individually, instead of targeting the entire country. This 
way each frequency can be reused in non-ad jacent areas and thus 
increase the system's capacity. The projected antenna would allow 
for individual control of the phase and power level of each array 
element. The objective is to produce a scanning spot beam that has 
low sidelobe levels, sufficient gain, minimum cross polarization and 
with minimum cross talk between beams. To aide in this design work, 
a considerable amount of computer modeling is needed in order to 
evaluate the effects of mutual coupling between the array elements 
on the above characteristics. This report presents a model developed 
to analyze the effects of mutual coupling in arrays of rectangular 
waveguides. The computer codes that were developed can be used to 
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aid in the design and development of modest array systems by solving 
for the radiation patterns, reflection coefficient and inter-element 
coupling coefficients. 

One of the problems that arises in the design of array antennas 
is that as radiating elements are placed close together the mutual 
coupling will change the radiation characteristic of each element. 

The effects of this coupling on an array are threefold: 1) changes 

in the total radiation pattern, 2) changes in the radiation impedance 
of the elements and, 3) variations in the polarization characteris- 
tics (_Amitay, Galindo and Wu, 1972). The pattern affects can be very 
pronounced. Sharp spurious nulls have been found experimentally in 
large arrays (ibid) at certain scan angles. 

There are several techniques or approaches that may be used for 
solving for mutual coupling and each has it's own strengths and weak- 
nesses. In this report the vector potential approach is used. This 
method is described in Chap. 2 along with a brief overview of some 
of the other approaches that may also be used. The vector potential 
method, as presented, uses the method of moments solution with vector 
potential modal functions as both basis and testing functions. These 
modal functions are used to approximate the equivalent surface magne- 
tic currents in the apertures. Once these currents are known, the 
radiation patterns and scattering parameters can be solved for. 

The numerical results from computer analysis of various sized 
antenna arrays with mutual coupling are presented in later chapters. 
The single aperture antenna is used to determine the number of modes 



3 


necessary for reasonable accuracy. The two and three aperture array 
results illustrate the coupling effects and are used as a basis for a 
first order coupling theory. Modest sized arrays are briefly looked 
at to consider the problem of coupling in larger arrays. Finally, 
results from experimental tests that were performed on a 2 element 
array are compared against computed data to verify the accuracy of the 
analytical work. Much of the computer data presented in this report 
deals with the testing of the first order coupling hypothesis. This 
hypothesis assumes that for a given array the mutual coupling can be 
solved for by assuming only coupling from a driven guide to it's un- 
driven neighboring guides and not back again. With this approximation 
it should be possible to reduce the computational time needed to ana- 
lyze large arrays. These results are presented here and will be used 
at a later time to formulate a technique for use in analyzing large 
arrays . 

The numerical and experimental results that have been obtained 
here are compared against both published and experimental results 
obtained by others. The validity of the analysis and the numerical 
results was thereby established. It appears that acceptable accuracy 
can be obtained using relatively few modes for the expansion of the 
magnetic current in each aperture. 

1 . 2 Review of Past Work 

A number of authors have studied mutual coupling effects in finite 
arrays composed of circular or rectangular waveguide apertures in 
a large ground plane. Mailloux (Mailloux, Jan. 1969) studied the 
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coupling between two rectangular waveguides using a vector potential 
approach. He assumed that no cross polarized aperture fields were 
excited. The self-admittance of a single waveguide and the coupled 
power between two waveguides were solved for and for some cases were 
compared against experimental data- 

The Fourier transform method (solving for the fields in the trans- 
form domain) was used by Borgiotti (Borgiotti, May, 1968) and then 
expanded on by Bailey (Bailey, 1972). The former obtained an expres- 
sion for the mutual admittance between two identical apertures in the 
form of a Fourier transform of a function related to the radiation 
pattern of each element. Bailey started with the same basic coupling 
equation as used by Borgiotti and developed a general expression that 
is applicable to apertures that are not identical in shape or excitation. 
He also included the case where a planar stratified region outside 
the aperture existed (i.e. a plasma layer). He computed the TE^ to 
TE^^ mode coupling between two circular apertures as a function of 
spacing. The reflection coefficient for a central and edge element in 
an array of 183 elements was also evaluated as a function of scan 
angle. Higher order mode coupling was treated in a limited way. 

Cha and Hsiao considered the coupling between rectangular aper- 
tures in terms of an aperture field expansion in TE and TM modes 

r r nm nm 

(Cha and Hsiao, 1974). The reflection coefficient for a central 
element in a 19 x 19 array was evaluated. 

Luzwick and Harrington also considered the problem of coupling 
between rectangular guides and used an expansion of the aperture field 
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in terms of TE and TM modes (Luzwick and Harrington, 1982). They 
present some numerical results based on a single basis function in each 
aperture with the aperture and waveguide of different sizes. The 
coupling between 2 identical apertures and the coupled power in a 49 
element array with a single driven element were evaluated. 

Steyskal (Steyskal, 1974) used a method similar to that of Luzwick 
and Harrington to solve for the coupling between circular apertures. 

The coupling coefficients for the center element of a 127 element 
array were evaluated and compared against measured values. 

Fenn, Thiele and Munk studied the coupling in rectangular aper- 
ture arrays using pulse and peice-wise continuous sinusoidal localized 
basis functions (Fenn, Thiele and Munk, 1982). Cross polarized basis 
functions were also included. The numerical computations were limited 
to the calculation of the input reflection coefficient for a single 
waveguide in an infinite ground plane, and for the edge and center 
elements of an 11 x 11 array as a function of scan angle for both E 
and H plane scans. 

Mailloux (Mailloux, Nov. 1969) studied the coupling of cross- 
polarized radiation between two square apertures and found that for 
certain orientations this is as large as the co-polarized coupling. 

The results for two aperture coupling were compared against experi- 
mental data. 

The above investigations have provided an excellent theoretical 
foundation for the analysis of mutual coupling effects. However, the 
numerical results were limited in both accuracy and scope and did not 
provide sufficient data for the design of high performance arrays. 


6 


Some of the major deficiences in much of this work are: 

1. Limited number of modes in aperture field expansions. 

2. Cross-polarization aperture fields are neglected. 

3. Lack of analysis of the beam forming properties of the array 
and the beam patterns. 

4. Lack of analysis of the effect of mutual coupling on the 
co -polarization and the cross-polarization of the beam. 

5. No consideration of the effects of mismatched waveguide sources 
on the radiation patterns. 
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CHAPTER 2 


ANALYSIS OF MUTUAL COUPLING 

2.1 Introduction 

This chapter will present an overview of various methods that have 
been used to solve for mutual coupling in antenna arrays that have a 
finite number of apertures. The method chosen for use in the computer 
model described here, the vector potential approach, will be discussed 
in detail (Sec. 2.4) along with a brief look at three other techniques 
that were considered. The overview of these three methods, included 
here for completeness, will include a short discussion of their respec- 
tive strengths and weaknesses, and the reasons why each wasn't chosen. 
The three are: Fourier transform approach (Sec. 2.2)» expansion in 

terms of TE or TM modes and piece-wise sinusoidal patches (Sec. 2.3) 
The latter two use the method of moments, much like the vector poten- 
tial approach, to solve the H field equations. The Fourier transform 
approach solves for the modal amplitudes in the transform domain. 

Following the discussion of the vector potential approach there 
will be a section on the needed mathematical techniques used to reduce 
the integral equations , along with one on the method used to derive the 
scattering parameters and radiation patterns. There will also be a 
section on a possible technique for approximating the mutual coupling. 
It solves for the first order coupling between a driven element and 
it's neighboring elements. Its use could reduce the computational 
time necessary to solve a coupling problem. 

Figure 2.1 shows the basic problem to be analyzed in the following 
sections along with it's coordinate system. A perfectly conducting 
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ground plane that covers the entire z=0 plane has a finite set of 
apertures cut in it. The apertures have H-plane lengths of a and E- 
plane widths of b. These apertures are waveguides terminating on the 
ground plane. The waveguides are individually driven by matched 
sources in the z<0 region. The z>0 space is assumed to be free of re- 
flecting obstructions. 

The assumption of matched sources will be relaxed later on. Also 
in practice a finite ground plane is used. The affects associated with 
a finite ground plane are not treated in this report. 


z 



Figure 2.1 General Problem 
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2.2 Fourier Transform Approach 

The Fourier transform approach solves for the aperture fields in 
the transform domain. The reason for using this method is that it 
results in some simplification of the equations to a form that is 
easier to work with from an analytical point of view. Since there are 
several approaches to the problem, a brief overview of one that appear- 
ed in the literature will be presented along with a general procedure 
that matches free space fields to the guide fields. 

The first approach is the one used by Bailey to solve for the 
coupling between circular apertures (Bailey, 1972). It starts by 
representing the transverse aperture fields as the sum of waveguide 
modal fields: 


E 

-t 


= y v * e ' + y v " 

•— • a — q q 

q q 


M 

e 


-q 


Bt-Z 


I' h' 

q ~q 


+ Y i" h" 

q q 


( 2 . 1 ) 


The h * and e ' are the TE modal functions and the h" and e" are 

-q -q ~q ~q 

the TM functions. Using the fact that the energy for each mode pro- 
pagates independently the equations are set up by assuming that each 
modal field is fed by a seperate source. This corresponds to an 
equivalent microwave network of N ports where N equals the product 
of the number of modes per aperture and the number of apertures. In 
matrix form the coupling is represented as : 

^NXl = ^NXN^NXl 


( 2 . 2 ) 


10 


from which the scattering matrix is readily found to be given by: 


[S] = [[Y q ] - [Yj] J[Y 0 ] + [Y]] 1 


(2.3) 


where 

IYq] = diagonal matrix of waveguide mode characteristic 
admittances . 

[Y] = coupling admittance matrix. 

I V ] = equivalent voltage vector. 

[I] = equivalent current vector. 


Once the admittance matrix [Y] has been determined, the scattering 
matrix can be found. 

The elements of the admittance matrix are found by solving the 
electromagnetic reaction integral (Harrington, 1961). 


i,3 - vir/Zd 1 1 • h d h 

1 2 S. 

l 


(2.4) 


where 


Y. , 


V. 

3 


= coupling between the i and j modes 


. th 


= surface of the aperture in i waveguide 


When the Fourier transform of the terms in the integral are intro- 
duced it is found that: 


Y 


1 


V.V. 
i 3 


MS 


[E 1 (k ,k ) X H j (-k T k )] 
— x y — x y 


a dk dk (2.5) 
_ z x y 


where E (k , k ) and H(-k ,-k ) are the Fourier transforms of the E 
— x y — x y — 

and _H fields. By solving for H and .E in terms of their vector 
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potentials, a double integral of a term composed of only the trans- 
formed E field results. After further manipulation the final 
integral for . must generally then be evaluated numerically. 

After numerically solving these equations for the admittance 
matrix, Bailey substituted it into (2.3) to find the scattering matrix 
[ S ] - Then using the S_^ terms (ifj), the coupled power between the 
array elements can be found. 

The second method illustrating a general approach to the mutual 
coupling problem using Fourier transforms is outlined below. Some 
of the steps and procedures outlined here can also be found in the 
references (Collin and Zucker, 1969). 

Initially, we consider a single aperture and define the trans- 
verse fields of the aperture in terms of modal functions in the fol- 
lowing form (where R is the reflection coefficient) : 

,nm (2.6) 


E (x,y)=V + in (l+R)e ._+ V " V“ e + T V, e, 
—a e , 10 — e , 10 e , nm-e , nm 2— h , nm — h , i 

- ™ n,m 


n,m 


where V , ,, is the amplitude of the incident dominant mode, and 
e, 10 r 

V , V, are amplitudes of the reflected TE and TM modes, also 

e,nm h, nm ^ 

E l | 

signifies the sum over all nm except n=l, m=0. The free space 
nm 

(z>0) field is a solution of: 


(V 2 + k 2 ) E = 0 


(2.7) 


which in the transform domain is: 

>2 


<k o - k t + 


az 


y) E(k »k ,z ) 
z — x y 


= 0 
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2 2 2 

where k = k +k and E (k ,k ,z) is the transform of E (x.y.z). 
txy — x y — ’ ’ 

This gives _E in the form: 


E(k , k ,z) = F(k , k ) e jk z Z > z * 0 
— x y — x y 


( 2 . 8 ) 


— 2 — T 

k -k • At z = 0 the transverse part F (k , k ) equals 
z 0 t —txy 

the Fourier transform of the aperture field expressed by (2.6). The 
z component of _F is obtained by using V*_E = 0 in the transform domain 
Thus , 

k-F(k ,k ) = k • F +kF =0 
x y ~t — t zz 


or 


-k • F (2.9) 

F = - A — 
z k z 

When we solve for the transformed H field using V x ji = we 

obtain 


-J k o z oS (k x ,y ■ -i 


( 2 . 10 ) 

After further manipulation the transverse part is found to be 


— a(k »k ) = 


-Y [k - k k • ] F x a 
0 0 — t~t — t — z 

k k 
0 z 


( 2 . 11 ) 


where F = E (k ,k ). The transverse magnetic field on the z = 0 plane 
— t —a x’ y 

is given by 

jk y y dk dk (2.12) 

x y 

The transverse magnetic field in the aperture may also be expressed 
in terms of the waveguide modes. When this is done and the Fourier 
transform is taken it is readily found that: 
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S (k ) = f { sin O* 1 ^) 
my l b 


W - - f I “ s o I 


C (k ) - -- j cos (~^) 
my I b 


(2.13) 


The /| J signifies the Fourier transform operator. The continuity of 
the tangential H field in the aperture then means that (2.13) equals the 
transform of H t (x,y) over the aperture, i.e. 


t if rr ej<kak,t ' )x +j 

^ j j j n j n 


^ k y k y^ y dxdy dk' dk' 
x y 


^tt J 0 J 0 

The function H (k ,k ) is given by (2.11) where F is the Fourier trans- 
— a x y — t 

form of E , which is given by (2.13) after setting all admittance para- 
— a 

meters equal to unity and taking the cross product with - a . The 
integration over x and y is readily carried out and leaves a two dimen- 
sional convolution integral still to be done. 

The above yields a system of equations with N unknowns (N = number 
of modes). The amplitude V“ determines the dominant mode reflection 
coefficient. In solving for these unknowns the aperture field can be 
viewed as a band-limited function of extent a along x and b along y for 
a rectangular aperture. As such the aperture field can be recovered in 


terras of the sample values at points spaced by ir/a along k^ and by 
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r/b along k in Che transform domain. We may assume that the aperture 
field is adequately represented by a finite number of modes and enforce 
the continuity of H at the wave numbers corresponding to k = nr/ , 

3 X 3 

= mr/b in k^-k^ space. A total of N spectral points are thus matched. 
This point matching procedure is equivalent to assuming that the free 
space field can be adequately represented by a spectrum of plane waves 
having transverse wave numbers that closely match those of the finite 
number of waveguide modes that are included in the aperture field ex- 
pansion. 

For the multi-aperture case it is only necessary to add to the 


free space term of each aperture, the sum of the Fourier transforms of 

the other aperture fields, each multiplied by e^x^x ~^y D y. where 

Da + D a is the seperation vector between the apertures involved, 
x— x y-y 

This concludes the brief overview of two Fourier transform ap- 
proaches, where the equations for the aperture fields or mutual admit- 
tances are solved for in the transform domain. The main disadvantage 
of the Fourier transform method is that the coupling integral (2.5) 
giving Y„ or the matching condition (2.13) involves the infinite k^-k^ 
domain. Thus the integral must be truncated or the matching must be 


limited to a finite part of the k^-k^ plane. In view of this, there 
does not appear to be any numerical advantage in using the Fourier 


transform approach even though the analytical results do have a certain 


degree of formal simplicity. 
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2.3 Expansion in Terms of TM and TE Modes and Sinusoidal Patches 
The approaches discussed in this section use the method of 
moments to solve the integral equation found by using the continuity 
of the H field in the aperture. Luzwick and Harrington (1982) use 
a single modal function per aperture as their expansion function. 

Fenn et al (1982) use piece-wise sinusoidal-uniform surface patches 
to solve a similar set of equations. A quick overview of these two 
approaches will be given here. Much of this material is taken 
directly from these two articles and the interested reader is refer- 
red to them for more detail. The equations and how they were derived 
are similar for both and start by looking at a single rectangular 
waveguide fed aperture in a ground plane. It should be noted at 
this point that both Luzwick and Harrington, and Fenn et al use a 
notation that is slightly different from that presented here. They both 

used M for J , F for A and Harrington used H imp in place of H *^' G . 

~ ~ m — ~ m — — 

First the equivalence principle is used (Harrington, 1961, Sec. 

3-5) to split the problem into two seperate regions by covering the 
aperture with an electric conductor. To keep the tangential compo- 
nent of the E field continuous across the aperture, equivalent magne- 
tic currents are placed on both sides of this conductor over the 
aperture areas as shown in Fig. 2.2. Since = H x j± , the currents 
on the half-space region are negatively directed. Next, due to con- 
tinuity of the tangential H the magnetic fields on either side are 
set equal. The waveguide region H fields in z<0 are 

H WG =H INC + H WG (J ) 

— t — t — t — m 


(2.14) 
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where 


WG 


' ( J ^ ) = tangential H due to the equivalent magnetic 


currents 

= tangential H due to the sources (J^ and _J^ ) under 
short circuit conditions at the aperture. 


The half space region then has: 


h hs = h hs (-j ) = -K HS (J ) 

— t — t m — t --m 


(2.15) 


Setting (2.14) equal to (2.15) at z = 0 yields: 


H WG (J ) + H HS (J ) = -H INC at z 
— t — m — t — m — t 


(2.16) 


When N-l apertures are added, we get for each aperture 


nf(j ) + T h“ s (j ) = i- 1 , 2 , ...n (2.17) 

— t — m . ' — t . — m . — t . 

1 j-1 1 3 


,HS, 


.INC 
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The H (J ) term is the radiated field at the i aperture due to the 

” C i ““J 


.th 


j magnetic current and is derived from the vector poten- 


tial A and the scalar potential 6 .. where: 
— m. . in 
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i._ // D .£WJ ds 

ii 4 irn J J j j r-r 
J 0 apert. J 1 1 


where 5 . = 


— ’ 


3 -ju 


H^ S (j ) = -j <o 
t . ti . m . . 

1 3 3-3 


VO 


(2.18) 



Closin 
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WG , 

In solving for the waveguide term H_ (J ) due to the J the 

— t — m — m 

two articles take different approaches. Harrington and Luzwick ex- 
press the field in modal form and solve for the mode amplitudes. 
The analysis is somewhat involved since they include the case where 
the guide dimensions are larger than the aperture dimension. Fenn 
et al use image theory to derive the field in their situation. Since 
both are quite lengthy, they will not be presented here. 

It is relevant, at this point, to present a quick review of the 
method of moments, since it is used both here and again in the vector 
potential approach. The method of moments allows the user to solve 
the following equation for f: 


L(f) = g (2.19) 

where L(), a linear operator, and g are known. The first step is to 
approximate f by: 

f = E » f (2.20) 

_ n a 
n=l 

where a = unknown constants 

n 

f = suitably chosen basis or expansion functions 

The f should be chosen so that if N= 30 the two terms in (2.20) are 
n 

equal. The substitution of (2.20) into (2.19) yields: 

N 

E a L(f ) = g (2.21) 

n=l n n 

Next a suitable set of N testing functions are chosen along with 
an inner product <a,b> . Taking the inner product of each W , also 
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called testing with , with both sides of (2.21) gives: 


a < W L(f )> = <W ,g>, m = 1, 2, . . . N 
n m, n m 

n=l 


N 

E 


^NxN^Nxl ^Nxl 


( 2 . 22 ) 


< V 1 L(f 1 )> < W x L(f 2 )> . . . 

< W 2 L(f 2 )> < W 2 L(F 2 )> . . . 



°1 


<W g> 
x > 




<W g> 

[a] = 

z 

[g] = 



• 

m • 


• 


Provided that [ i ] is nonsingular then 

[a] = m _1 [g] (2.23) 


where 


M 


NxN 


Proper choice of f W and <a,b> dictates the accuracy of the solu- 
n, n 

tion. A more complete discussion of the method of moments is given 
in Harrington's book (Harrington, 1968). 

The basic equation (2.16) is in a form that is readily solvable 
by the method of moments for the unknown (which is related to 

the aperture field JS) . Here again the two groups differ in the way 
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c’nev approximate J . Harrington and Luzwick use full domain functions 
analogous to the TE and TM modal functions with the IE mode functions 
having the form: 


_ n ,mrx, . ,nTTv. m . / m~x N ,nmv,-| 

C fa -r- cos(— — ) sm a — sm (— — ) cos (— r^)J 

L — x b a b — y a a b -* 


In the work presented in their article only a single function was used. 
Fenn, et al used piecewise sinusoidal-uniform basis functions of the 


form: 


a K 
— x n 


sin [3 ( n L- I x I ) ] 
sin(6I.) 


(n-l)L-'x< (n+l)L,0<v<W 


for overlapping surface patches of length 2L and width W as shown in 
Fig. 2.3, where 3 = 2i/\ and n = basis function number. 



Figure 2.3 Piecewise Sinusoidal Basis Functions. 
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These can be used for both x and y directed currents. Seven sinusoids 
per aperture for rectangular apertures was found to give satisfactory 
results . 

For the choice of testing function both Luzwick and Harrington 
and Fenn, et al used the same function as they used for their basis 
function. This technique is known as Galerkin's method and in the 
former case would allow for the use of the orthogonality of the modal 
functions to reduce the complexity of the integral equations. For 
the inner product both groups used: 

<a,b> = ff dS 

where the integration is over the aperture. 

Once the II] matrix .is solved for numerically it can be inverted 
and the amplitudes of the expansion functions determined. With these 
amplitudes the approximation for the equivalent current is complete 
and the scattering parameters can be found. 

Fenn et al used their approximation to solve for the reflection 
coefficient of driven elements in several array situations while 
Harrington and Luzwick found the coupled power (20 log I S ^ I ) between 
apertures in several configurations. Using the piecewise sinusoidal 
expansion approach, Fenn et al were also able to look at the edge 
singularity present in the aperture fields. Some of these results 
will be presented in Chapter 3 as a comparison with the results found 
with the method used in this report. 

A good general reference covering the basic ideas presented here 
is the report by Harrington and Mautz (Harrington and Mautz, 1976) 
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It should be noted that the definition of the mutual admittance 
parameters is quite arbitrary. It can be chosen as the coupling co- 
efficient between the n C ^ basis function in the aperture i and the m C ^ 
testing function in the aperture j. The various choices of basis and 
testing functions used will then give different expressions for the 
admittance parameters. As a general rule, the scattering matrix re- 
lating the amplitudes of the propagating modes in the waveguides is 
a desired end result so the evaluation of the admittance parameters 
is often only an intermediate step. 

The main disadvantage of the TE and TM mode technique is the fact 

that each coefficient will have to be split into an a and an a term. 

—x — y 

As we will see in the next section, the vector potential approach 
overcomes this difficulty. The reason for not choosing the local 
basis functions is that the size of the matrix that would have to be 

inverted for arrays of square apertures would become very large for 
modest sized arrays. This would severly limit the cases that could 
be solved. 

2.4 Vector Potential Approach 

The vector potential approach is very similar to the methods 
discussed in the last section. The main difference is in the modal 
expansions used to determine the fields and the approximation for the 
magnetic currents. The approximation uses vector potential modes for 
basis functions. This allows the equations to be solved by first find- 
ing the magnetic vector potentials, which are more straight-forward and 
easier to work with. The TE and TM modes in the waveguides can be 
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found later if this information is required. 

Consider again the single aperture case (see Fig. 2.2). There 

exists a J (magnetic current) on both sides of the aperture surface 
— m 

equal to a x _E , where _E is the aperture E field. The form of this 
Z 3 a 

current using the vector potential modal functions (also called a hy- 
brid mode set; Harrington, 1961, Sec. 4.4) is: 

n;[ m^ n 2 m 2 


J 

— m 


L L 


nm 


„„ 4 EE 0 a 

n=0 m=l nm 


nm 


(2.24) 


where 


/On^m. 1 ^ sin (— ) 
^nm = ( 'T^ ) 


COS 


( m -f) 


a 

—x 


(2.25) 


* / £ 0n £ 0m x ^^2 ^-rry^ ,mrx s a 

9 = ( r — ) sm (— 7f-) cos ( ) — y 

— nm ab b a 


These modal functions are orthogonal so that: 


ff* ip , , 

J I nm n m 


ds = 0 


for n^n' or m^m' 


This relation also holds for <J> and for any combination of the two 

nm y 

sets of functions. 

The next step is to set up the integral equations using the 
magnetic vector potential. _E and H are found using 


e_E = - V x A 
0— — — m 


J"Vo - ” < k o + ™->4. 


(2.26) 
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The continuity of H across the aperture gives: 



+ _2Y* ) ( A 


WG 


+ a inc ) = 

“Til 


(V 2 + 7 7 * ) 



(2.27) 


where 


,WG 

A = 

— m 

a hs = 
— m 

( INC 
A = 

— m 


Vector potential due to inside guide. 

Vector potential due to - J radiating into 
free space. 

Vector potential due to the sources under 
short circuit aperture conditions. 


These three vector potentials are the quantities that need to be 

solved for in terms of the J . 

— m 

Using the magnetic vector potential integral (Harrington, 1961 
Sec. 2.10), duality, and image theory: 



(2.28) 


where R = |._r-_r 1 | (r; = field vector and _r ' = source vector). 

For the vector potential in the waveguide due to the it's Green's 

function needs to be found. Using duality again the form of the 

\ 

equation to be solved is (Collin, 1960 Chap. 2) 


which gives 

Since there are 
function should 


<v 2 * 


(2.29) 


«• ■ ffi 

o* / o 


J dv ' dx ' 
— m 


only x and v directed currents the required Green's 
be a dyadic of the form: 
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G = G a a + G a a (2.30) 

= xx x x yy y y 

where (V 2 +k 2 ) G = -2 £q 6 (x-x') 6 (y-y 1 ) 6 (z ) (2.31) 

(V 2 + k 2 ) G yy = -2 e 0 6(x-x')S(y-y')S(z ) 

For G we may assume an expansion of the form 


G 

xx 


£ 

n, m 


G ip 

nm nm 


-jS 

e nm 


(2.32) 


where the G ^ are unknown amplitudes (Collin, 1960, Chapter 2). 
When we substitute this into (2.31) we obtain 


£ 

n,m 


G * [kj * ( I T L ) 2 ’ O 2 

nm nm U a b 


= - 2Eq <5 (x-x') 6 (y-y') 6 



(z) 




Iz ! 


nm 


(2.33) 


Multiplying this by A (x,y ) and integrating with respect to x and 

nm 

y gives : 

2 

G (g 2 + |— 2 ) e ^nnJ 2 ^ = -2e_ <p (x',y') <$(z) (2.34) 

nm nm d z 0 nm 


Integrating this equation over z from 0“ to 0 + gives 


9 e 
nm 3z 


-jgnm| z | 


0+ 


0 - 


= -2e ip (xjy ' ) = -2 j 3 G (2.35) 

0 nm nm nm 
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so finally 


G xx ^ 77 — ti (x,y) <1> (x,y ) ej3 nm 

n,m 16 nm nm 

’ J nm 


(2.36) 


Similarly G is found Co be given by: 
w 


’v v = y' — — <J (x,v) <{> ( x y ' ) e J °nm' 
-jg nm nm 

n , m J nm 


. i „ t \ jS 


WG 


(2.37) 


Therefore, Che A' ac z = 0 is: 
— m 

■a -b 


(G C $ (x jy ' ) a + G D 4> (x',v’) 

xx nm nm } — x yy nm nm ’■ 


a wg = l f r 

m n,m *^0 *70 

a ) dv’dx' = £ (x,y) + £ D * (x,y) (2.38) 

— v • ^ id — nm *■“' i S nm nm 

1 n , m J nm 


n,m nm 


INC 

Finally, for Che A ‘ under shore circuic condicions we have 


.INC . , , -jS lA z , j3 1A z. 

A * A in 'i in (e J 10 + e J 10 ) 

— m 10 —10 


(2.39) 


where A^q is a known amplicude consCanC. 

INC 

This makes E_ a TE^q m ode chat is zero at z = 0. 

The substitution of (2.39), (2.38) and (2.28) into (2.27) gives 


<k o l + ^t V/* f -r [ = +£ o 

J o J o o 


a + a a ) ] 
c-x -y— y 


d^y^dy'd*' . CkJ I 4- Vt )-(-2A lollo ) ^ 


(2.40) 


where the V is used since only the tangencial fields are being 
analyzed (z has been set to zero), and I is the unit dyadic. 
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This is the integral equation that is solved by the method of 

moments. A finite expansion of J in terms of and $ is an 

— m — nm — nm 

appropriate approximation for use in this method. The basis and test- 
ing functions are chosen as the same ip and $ . The following 

— nm — nm 

inner product is used: 


<a,b> 



dS - 


(2.41) 


where the integration is over the aperture opening. 

This choice of inner product, basis functions and testing func- 
tions yields a set of admittance parameters for the elements of the 
[ £] matrix. This can be seen by looking at the integral equation 
that results for these elements. 

<J 2 ,H L > - f f J 2 ■ H 1 dS 
-m - JJ ~m ~ 

S 

a 


where 


= H field in aperture //I 

2 

J = Current of aperture #2 
— m 

= Area of aperture It 2 


By definition, this is the reaction integral. (Harrington, 1961, 

Sec. 7-7). It is related to an admittance by (2.4). So after testing 
(2.40) with the j and £ the resulting matrix equation will be of 


the following form: 


[Y][V] = [I] 


with the proper choice of units for C and D (volts) . 

nm nm 
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HQ WG 

The [Y] is, in fact, the sum of two matrices, [Y “ ] + [Y ]; 

HS \<G 

each due to the respective magnetic field (H and H ) that is being 
tested. The [I] is the tested incident magnetic fields, while the 
[V] is the amplitude matrix that needs to be determined. 

In order to make the documentation and further manipulation of 

the tested eqns. easier, the following admittance parameters are 
introduced (where 1 and * are the testing functions): 


pq ,nm 


Coefficient of C testing <i> against 'ii 

nm ~pq ° — nm 


pq ,nm 


Coefficient of D testing y against J 

nm — pq ° — nm 


pq ,nm 


Coefficient of C testing $ against 

nm — pq ° — nm 


Y y 1 

pq ,nm 


Coefficient of D testing i against : 

nm — pq ° ~r r i 


Using these admittances in writing the resulting equations we get: 


i r l 

p , q L a » m 


(Y X ' X 1 C ) + H (Y X l y l D ) « :L ^ L0 
pq,nm -nm n,m pq,nm nm 


(k 2 - 11 ) 5 5 
K 0 a 2 ; pi qO 


^ (Y y : 1 C ) + Z- (Y ■ - - D ) = 0 (2.42) 

, m pq,nm nm n,m pq,nm nm 
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where jLI is over the 4 values and is over the ® values, 

p.q , p q p»q pq 

and likewise for / - and JLj • 
n , m n ,m 


A typical term in the tested equation is of the form: 

f * 2.„ 2 • A (J ) dS (2.43) 

J J nm t t — m — m 

HS 

where is given by an integral in certain cases ( A^ ) . To make 

these integrals easier to carry out, we use the following integration 
by parts procedure. Consider: 


7 ’ , b V -A ) = (V * )(7 *A ) + \p -V 7 • A 
— t ' -^nm— t m' — t -Tim — t — m -Tim— t— t — i 


— m 


(2.44) 


which when applied to (2.42) yields 


/ n. w V -A dl - ( ( 7 • ip 7 -A dS - 

— — nm — t — m J J — t ^im t m 

C 

- /i (*■ A )7 -ib dl + // 

J m -t -nm JJ 


ib dl + / / (V V -ib ) -A dS 
— nm / / t - t — nm — m 


(2.45) 


Note that n «A W ^ = n*A^^ = 0 on the boundary. 

— m m 

Using (2.42), (2.45) and the orthogonality of the modal functions 

in solving for the four admittance parameters we get: 

2 p 2 tt 2 


- i(k“ - £ — 7 — ) o 2 

Y X 1 X 1 = _ — Q SZ 5 6 +f k 2 _£ 


pq ,nm 


pn qm 0 a~ 


pq 


* ff 

J o J o 


■b (:: 

D q 


■" [ Jff 


e-j k 0 R 


b (.x 1 , y') dv'dx' 


2ttR nm 


^ dy dx ^ 


3 ib 

.. m (x , 

o X 
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x 1>' l 


a =' >. ■ - 

r n V-> [ /? 

Jo Jo L •'o J o 


= -UE3I1 0 0 - ££ 


pq,nm a'oS pn qm ab 

P<1 


-jk n R 
e J 0 


2~R am 


*(*’./) dy' 


y'dx' ^ dy dx - J>" hr 

[ f 0 £ 4 - dy ' dx '] 


(x,y) 


(2.46) 


(2.47) 


Y yi x i = Y Vl 


pq ,nrn 


pq ,nm 


with che <i/ and O's interchanged and Che — 

c X 


replaced bv 

'ey 


( 2 . 48 ) 


Y y iy i 

pq ,ntn 


f x l x l 
pq ,nm 


with the * andi>'s interchanged, the h— 

cX 

g 2 n 2 - 2 

replaced by -r^j and the two (k^ - • ~ 0 's replaced 
^ a 2 * 2 ’ . 

by (k n - ^"2 ) s (2.49) 


The next step in the analysis is the multi- aperture problem. The 

two aperture case is shown in Fig. 2.4. Let C 2 and D 2 be the 

nm nm 

modal amplitudes of aperture "1 and C ^ and D 2 be the amplitudes 
r run nm 

of aperture «2. The currents for aperture *2 are then: 


J 2 * 
— m 


n 1 mi 

E E 

n=i m=u 


n? m 2 


C *" ZD 2 * (X ,>’ ) 

n m “nm l l < nm — nm 2 2 

n=u m=i 
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INC 2 

Also the A of •'•2 is: 2A ■!) (x ,v~) at z = 0. If we define 

— m 10 —10 2 2 

Y as the coefficient of C 1 for b (x-,v ) tested against 

pq , nm nm nn / 7 * / ° 


nm 


* \ -‘•o » ' O ' 

pq 2 ’'2 


b (xi .vi ) and the remaining aperture to aperture admittances in a 
nm 


similar manner the resulting tested equations are: 


■ ■ 

W2 

E 

n, m 

xx , xx 9 , xy , 

(Y 11 C H 12 C 2 ) + V (Y 11 D 1 + 

nq , nm nm pq.nm nm i—u pq.nm nm 

n,m 

x y 

Y 1 2 
pq , nm 

D 2 ) 
nm 

- 2Al0 £ 5 (V 2 n 2 / 

" d pl 5 q0 (k o -^ ) 


r 

p,q 

[ z 

y x . y x ~ 

(Y 11 C 1 + Y 1 2 C 2 
pq,nm nm pq,nm nm 

y y i 

) + y (y 1 1 d 1 + 

L~t pq.nm nm 

n,m 

y ' y 2 
Y i 2 
pq ,nm 

D 2 ) 
nm 



Z 

p,q 

f E 

L n > m 

X X . 1 X 0 X- 9 

(Y 1 C " + Y C “ ) 

pq,nm nm pq,nm nm 

' x 2- v 

+ Z— < ( Y 3 - 

n,m pq , nm nm 

X V 

Y 2 ' : 

pq ,nra 

2 a 2 ) 

nm 

~ 2A i0 2 , , 2 n 2 

- E 0 5 pl 6 ,0 <k 0 - a 7 >J 


E' 

p,q. 

’ E 

n ,m 

y x i y x 9 

(Y 21 C 1 + Y 22 C 2 ) 
pq,nm nm pq,nm nm 

+ y (Y 2 1 D 1 + 

pq.nm nm 
n , m n 

y y 9 

Y 22 D 2 

pq,nm nm 


(2.50) 


where Y 


x.x. 

l l 

pq,nm ’ 


x .y . 

{ 1 1 
pq ,nm 


etc. do not depend on whether i 


1 or 2. 




a 


* n 


Ill- 



Figure 2.4 Two Aperture Problem. 


Using R = V -D x ) 2 + ( yi -y' 2 -D y 


and solving for the new 


admittance parameters we find that: 


pq ,nm 



/, 2 p 2 t 2 , , 

(k 0 ' } 




2 rrR 


'4__ (x ? ' y ' ) dy\ dx ? ' 
nm I , i • *_ z 


dv 




// 


l) (x n y ) 
DC 1 1 


0 X 


u; 
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omeiNAL P Aag is 


-jk R 
e J 0 

2ttR 


• - nm , y • ) dv ’ dx' 
— nm 2 9 2 2 


dl. 


(2.51) 


f X i y 2 = -pq - 2 

pq ,nm ab 


^/T 


2 2 


? pq < *l ,S 'l > 


rr 

J ojo 


/ c'b (x ,v ) f f 

J J 

1 . * / 0 J 0 


e _:ik O R 

- — $ ( x ' ,v' ) 

2ttR nm 2 ’ 2^ 

a r'o 


e - jk o R 

2ttR 


* - 4m (X 2 * y 2 } dy 2 dx 2 dl l 


(2.52) 


yi *2 

Y 

pq ,nm 


x lY 2 

Y with the ij;'s and i's interchanged and the 

pq ,nm ° 


— - replaced bv -r 2 - 

0 X 1 ’ 3 vj 


(2.53) 


Y 


y y A 

1 2 


pq ,nm 


Y with the ip' s and <5>'s interchanged, the 

pq.nm 

.. 9 2 tj 2 

replaced by — and the (k 2 - £-- 7 — ) replaced bv 
cx^ o U a 


( 



(2.54) 


For the other 4 terms reciprocity shows that: 


x.y . 

1 J = 
pq ,nm 

y.x. 

Y J 1 
nm,pq 

y.x. 

1 J = 
pq ,nm 

x.y. 

Y J 1 
nm, pq 
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X . X . 

X . X . 

■ 1 J = 

Y J 1 

pq , nm 

nm, pq 

y,y. 

y.v. 

' 1 3 = 

Y J 1 

pq ,nm 

nm, pq 


(2.53) 


By using(2.45) - (2.55) it would be a straight forward task to extend 

this analysis to the cases of three or more apertures. So those 
situations will not be presented here. 

These equations completely describe the system of several aper- 
tures in a ground plane and can be solved directly for the C 's and 

nm 

D ' s There is still one remaining problem with them and that is the 
nm 

three and four-fold integrations that need to be carried out numeri- 
cally. If Simpson's rule is used and n or m were allowed to be as 
large as 3 then the minimum number of intervals for an acceptable 
approximation would be 6 (4 points per period of sinusoid being in- 
tegrated). So, if 8 were used, for somewhat better accuracy, then 

4 

each four- fold integration needs 8 ~ 4000 operations, where each 
operation involves both exponentiation and trigonometric computations. 
Considering that each guide can have several modes and there are 
several guides and the number of terms to be integrated before sym- 
metry is considered is the square of the product of the number of 
modes times the number of apertures f this rapidly leads to an enor- 
mous amount of cou time. So the integrals need to be reduced as far 
as possible. A double integral using 20 points would only use 1/10 
the computing time, so by reducing these integrals to double integrals 
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fairly accurate results can be had at a fraction of the cpu time. 

The reduction of these terms analytically will be outlined in the 
next section. 

At this point some of the advantages of using the vector poten- 
tial method and an expansion of it can be readily appreciated. First 
of all there is a single set of amplitude coefficients along with the 
x and y components of the vector potential. In contrast an expan- 
sion of the aperture field in terms of TE and TM modes results in 
two sets of amplitudes associated with each component of the electric 
field. This results in a splitting of the admittance parameters into 
sets representing coupling between TE-TE modes, TM-TM modes, and 
TE-TM modes. As a result the notational complexity is increased sub- 
stantially. In the vector potential approach no distinction between 
TE and TM modes needs to be made. Luzwick and Harrington use both 
the vector and scalar potentials and do one integration by parts 
(Luzwick and Harrington, 1982). In the approach given above, two 
integrations by parts is carried out which leads to simpler final 
integrals to evaluate. 

2. 5 Reduction of the Three and Four-Fold Integrals 

By making a change of variables it is possible to reduce the 
three and four- fold integrals encountered at the end of the last 
section to double integrals. This greatly eases the computational 
costs by reducing run times and reducing the needed programming tasks. 
The first step will be to analyze the coupling terms between modes 
within a single guide. The steps used to solve this case can then be 


36 


applied to the coupling terms between different guides. 
Consider the following integral: 


r a f bfa f b 

I 

J c J 0 J 0 J 0 


g(x ,x ' y,y' ) <l> (x',y' ) dvdxdv’ dx' (2 

nm 


where : 


g(x,x,' , y,y') = 


££V 

2ttR 


R = [(x-x' ) 2 + (y-y’ ) 2 ] 1/2 


Let 


u = x-x 


and 


= x+x' . This gives x = 


v+u 


and x 1 = 


v-u 


and 


dxdx ' 


c x 3 x 


l/ 2 1/, 


3u 3v 




3x! 3x' 

dudv = 

-l/ 2 l/ 2 

dudv = — dudv 

3 u 3v 





The new limits of integration are: 


X 

il 

o 

corresponds 

to 

V = 

-u 

x = a 

corres ponds 

to 

V = 

2a- u 

o 

ii 

-x 

corresponds 

to 

V = 

u 

x f = a 

corresponds 

to 

V = 

u + : 


Thus the mapping into u-v space is a parallelogram as shown in 
Fig- 


. 56) 


2.5 



Figure 
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Now lee v-v = w and y+y' = 2 , so that y = 


w+z 


and v ' = ~ . 


2 •' 2 

A similar mapping cakes place here. These ewo changes, when made to 
(2.56) yield 


1 / / / / ,v+u w+z 


) • g ( u , w) (^ 7 -’ ^ 7 ^) dudv dwdz (2.57) 

nm 2 2 


Equation (2.57) can now be integrated analytically over v and z. 

2 2 . 1/2 

The g(u,w) is an even function of u and w since R = (u + w y 
while the other term in (2.54) is: 


( £ 0n''0m £ 0p~ 0q) 

2 . 2 

a b 


U2 


sin ( ■ ( u+v ) ) sin ( 7 - (v-u) ) cos (-f -r (w+z)) 

2 a 2 a 2 b 


cos (^(z-w)) 


(2> 58) 


To integrate out the v term, start by looking at: 


sin 0 ^(u+v)) sin (^(v-u) ) = sin (-J^pO cos Gy^ 1 ) sin ( 7 ^-) 


. D7T V , . ✓ yu, . /y v x ,111‘Vs / n Ll , 

cos (- 7 — ) - sm — ) sin (-= — ) cos (t — ) cos ( 7 — ) + cos (5; — ) 
2 a 2 a 2 a 2 a 2 a 2 a 


ITU, 


,nfu. 


.2TTV 


,nirv. 


,n^u, . /Pfrv. . ,mv. . ,n^u. ,pru>. .nTiv, 

cos (-y^ - ) sin (+ 7 ^-) sin ) - sm cos cos ( 'Ta 


sin (f=L) 
2 a 


(2.59) 


This is multiplied by g(u,w) and integrated. Since g(u,w) 
is even in u, only the terms that are also even in u will contribute. 
Expanding these terms gives: 

7- j [ cos (-^— ■ ttu ) -cos (^7^ ttu)| [cos(-^-^- TTV ) + COS irv)J 

4 2a 2a 2a 2a 
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+ [ cos ( 2 ~) 77 u) + cos ( 2 ” u u)| f-cos (' 2 ^' 1Tv ) + cos 71 v )j| 


(2.60) 


The integral over v is then over the strip shown in Fig. 2.6 


v 



Figure 2.6 Integration Limits 


Since it is an even function in 


in u let / | | du = 2 / 

J-a J 0 


J 


du . 


The only terms left to be integrated are cosines, where a typical term 
yields : 


-»2a-u 

L c 


, . . . , , .1 + COS (2as) . 

cos (sv)dv = - sm (su) f ] 


(2.61) 


where s = ir or it- Using (2.61) in integrating (2.60) over 

2a 2a 


v leaves: 


cos [(£ ir 7U] sin t l£ g } ,u ] [1 + cos ((p+nh)l 

(P-Kl) * 

2a 
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-cos 


[ ( -P^^ 7TU ]sin[ -^^ " u ] [ l-K:os((p-n)T) ] | T^Cu) a 


za 


(p— n) i 
2a 


V ~0n - 0p 


(2.62) 


2a+u 


Note that integrating over the LHP strip in the u-v space ( 


L 


n 


du) 


would give T^(-u) for u< 0 So the result in (2.62) should be 


V 


£ 0n £ 0p 


T 1 (|u|). 


For the cos (-^ (w + z)) cos (■— (z-w)) term the even functions 
Id Zb 


in w are : 


■J | COS ( * W) COS ( ~^j~ * 2) + COS (-^j^-1TW) COS (-^2t!Ill rj z ) 


(2.63) 


Integrating (2.63) over z using (2.61) yields: 


j J- cos tt w) sin (^LT^ t t I w i ) [ 1 + cos ((q+m)^) 


2b 


(q-Hn) tt 
2b 


cos 


( H w) sin ( Vk^ wl) [ 1 + cos ((q-m)^) 1) T 0 ( [ w | ) b 


2b 


(q-m) 

2b 


For the integral 


^ //// 


V £ 0m 'Oq 


(2.64) 

( H±Z). o(u w) $ (1±H 3HL) 

2 > o'-u.w; 2 , 2 > 


dudv dwdz the two equations that we get after integrating over v and 
z are : 

T^(juj) = ( | o | ) with b replaced by a. 

(2.65) 


t 4 (|wD = VI w | ) with a replaced by b. 
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Let 


Vr - V % L 


+ cos [— — ] [ 1 +cos ( (p+n) n) J s in 

(p+n) ir 
2h 


cos [ l£±^Il! l] s in [XE_S12 lLlI] ji + C os( (p-n)rn) ]| 


( 2 . 66 ) 


(p-n)ir 

2h 


Then 

T x ( |u|) 

= T + 
pn 

(u/a) 


t 2 (M) 

= T 

qm 

(w/b) 


t 3 (|u|) 

= T 

pn 

(u/a) 


T 4 ( | W l ) 

= T + 
qm 

(w/b) 


(. 2 . 67 ) 


Using the same techniques and variable changes the triple integrals 
can be reduced also. This then gives: 


? 7 

A 2 n^TT^ x 

Vi -MV" 

Y = 6 6+ (k_ 

pq.nm B pn qm 0 

nm 


' 2 TT 2 

’" Z ~ ' 2tt JO •'O 


2 p 2 tt 2 ^ 1 
- — ; — 


/ / 

Jo *'0 


+ _ e -jk Q (w 2 +u 2 ) /z 

T (u/a) T__ (w/b) — 7 7 jy — dwdu 


pn 


qm 


(w 2 +u 2 ) /2 


■f? 


(e 0n £ 0p ) 


2 


2a" 


L-j^Cx’^-fw 2 ) /2 e - j^ 0 ((a - x' 1 ) 2 +w 2 ) 1/2 

- cos ( pn) 


(x^ 2 + w 2 ) ^ 2 


9 1 / 

[(a-x^)+w j '2 
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sin • T (w/b) dwdxl 

a qm 1 


( 2 . 68 ) 
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Y 11 _ _21 JUL 6 6 _ 

pq,nm 6 ab pn qra 2ab 
nm r 
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Jo Jo 


T (u/a) T (w/b) 
pn qm 


e -jk 0 (» 2 +u 2 ) l, Z 


(w 2 +u 2 ) * 


dwdu + 


rf 

~'o -Jo 


'b (£ £ ) 1 /2 

v ^ Qm ; 

2ab 


2 / 
3 -jk.(u+y. 2 ) '2 


Oq Om o r e J 0 J 1 

TZJT 


(u 2 +y’ 1 Z ) V 2 


- cos (qir) 


-jk 0 (u 2 +(b-y i ' ) 2 ) 1 /2 


n ] sin C 32 ?*) T" (u/a) dy! du (2.69) 

2. 1/ b pn 1 


(u 2 +(b-y i ’)‘) */2 


f y l X l 

pq ,nm 


nm^r 


_ 22l 


3 ab ‘"'pn ^qm 2ab 
nm v 


ff 


T + (u/a) T ~ (w/b) 

pn qm 


e -jk 0 (u 2 4w 2 ) 1 /2 dudu + f a r b ( £ 0n !0 £) V 2.CL 


, 2 . 2 . l.„ 

(u +w ) /2 


IT 


0 •'O 2ab 


-jk n (x' 2 +W 2 ) L/ 2 
e J 0 i 

(x^ 2 +w 2 ) 1 >l 


, . -jk n ((a-x.) 2 + w 2 ) ^2 

cos(pTr) e J 0 1 _____ 

( (a-x|) 2 + w 2 ) 


sin (3221^1) . x (w/b) dwd t 
a qm x 


(2.70 ) 


,¥i -i p.V, 

pq ,nm b nm 0 ,2 
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5 6 + (k 2 - 

pn qra 0 ,2 


Sliof\n { 


(u/ a) 


/ 2 o, ^ X/ 9 
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. 2 ^ ,2.1/ 

-jk (u +y^ ) 2 cos (q 


(u 2 ^i 2 ) V 2 


2 2 / 2 

\ -jk.(u +(b-y 1 ' )~ , 

!f) e J 0 - 1 y . .miTyi . 

l sin (-—) 


(u 2 +(b-y’) 2 ) 1/ 2 


T (u/a) dy' du 
pn 1 


(2. 71 ) 


These are the Integral equations that are to be solved numerically. 

It should be noted that these four integral equations must 

he further manipulated before applying any of the standard numerical 

integration techniques. First, if p=n in (2.66) then there exists 

a singularity in one of the terms where l/(p-n) = ^/n. The T (t/h) 

fa Pn 

is incorrect at that point. The reason being that / cos (0 u) du 

JO 


should be 


/ 1 

Jo 


du = a not 


sin 


(Ou) 


'0 

So for these special 


0 


cases the p+n or p-n equal to zero should be substituted in first 
before integrating over the v or z . This gives: 


T 

0.0 

(t/h) 

= 2/h(h- |t !) 


+ o 
o 
H 

(t/h) 

= 0 



T + 

(t/h) 

1 

\ [cos(2piT) + 1] sin(D*ltl/h) 

(pTT t/h) ( 

PP 

h i 

i — ' 2 cos 

pir/h 

T “ 

(t/h) 

= i 1 

r.os ( 2h— h— sinf^Hfl + COS 

( 2 p TT )1 1 

PP 

h I 

. ' h ' pVh J 


With this change in (2.6o) a new problem arises. That of what to do 
about the singularities that now occur in the integral equations at 
u = w = 0 for the three non-zero cases listed above. The term is : 


44 


II 


r - (u) T ± <w) e-^0<“ 2 *' 2 > 1/2 
pn q m ' 


dudw 


(2.73) 


(u 2 +w 2 ) ^2 


and when integrated numerically the u = w = 0 point is a singular point 

since T — and T — are no longer equal to zero in these cases. To 
pn qm 

get around this difficulty, the singularity must be subtracted out. 

This is accomplished by adding and subtracting 


II 


T -(0) T— (0) dudw 
pn qm 


from (2.73) leaving: 


VT 

> u 


4* w 


ii 


T - ( u ) T - (w) - T - (0) T - (0) 
pn qm qm 

(U 2 + W 2 ) V 2 


dudw + 


If 1 “ 


(0) T7(0) 

qm dudw 


2 +w 2 ) ]/ 2 


The first term then equals zero at u = w = 0 and the second term's 
denominator can be integrated analytically to yield: 

b In (a + (a 2 +b 2 ) ^2 ) + a In (b + ( a 2 + b 2 )^2) -b In (b) - a In (a) 


Note that for p # n in (2.66) that the numerator is zero in the 
integrals given in (2.73), so no singularity exists. 

The last step then is to reduce the equations for coupling 
between different guides. First, it should be noted that (2.66) is 
no longer sufficient to cover these 4 cases since g(u, w) is 

2 

no longer symmetric in u and w, since a typical R here is ((u-D ) + 

2 V 2 + X 

(w - n ) ) . So a second term needs to be added to the T — (t/h) 

y pn 

to include the non-symmetric terms dropped earlier. The (2.59) non- 
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( 2 . 78 ) 
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X l y 2 
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-p^ 2 q f a 
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1 , P n 
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(u/ a^ t 

(w/ b ) + 
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4ab JL a J 

R * (w/b) ] 
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b x z y x z y 


[ T (u/a) 4- R (u/a)] dy' du (2.79) 
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y x 2 fa fb 
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G(u-D x ,w-D y ) dwdu + ^ ( £ 0n £ 0p )1/2 


l:f; 


nirx , 

a . , 2 \ 

sm ( ) 


( G(0-x'-D , w-D )- cos(pTr) G(a-x'-D , w-D )][ T (w/b) + R (w/b)] 
2 x y 2 x y qm qm 



R + (w/b)] G(u-D . w-d ) dwdu + 
qm x y 

[ G(u-D , 0-y ' - D ) - cos (qir) 
x z y 

R (u/a)] dv* du 
pn • 2 


(2.80) 


[ T (u/ a ) + R (u/ a )] [ I + (w/b) + 
pn pn L qm 



G (u-D , b-y ' -D ) ] [ T (u/a) + 
x z y pn 


(2.81) 


, 2 2. If. 

e -j k 0 (u +w ) 2 

where G(u,w) = ■ — -= — x — — 

2ir(u +w ) 


It should also be noted that in looking at the equations, it can be 
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seen that 


t Vl 

pq ,nm 


y l x 2 

I 

nm, pq 


(2.82) 


For aperture to aperture coupling there is no need to worry 
about singularities since R will never equal zero. The only way it 
can happen is if the distance between the apertures exactly equals the 
size of the apertures with zero wall thickness, but then there wouldn't 
be two apertures only one larger one. 

This concludes the mathematical manipulations needed to reduce 
the integral equations. The terms are now all double integrals that 
can be easily evaluated by Simpson's rule. 


2.6 Scattering Parameters and Radiation Patterns 

By using any of the techniques described in the earlier sections 

the amplitudes of the equivalent magnetic currents (or of the E in 

the case of the Fourier transform approaches) can be found. Once 

these are known the information must be put into a form usable by the 

antenna designer. The parameters that can be solved for using these 

quantities are the scattering parameters and the radiation patterns 

(which provide the side lobe levels, cross-polarization, etc.). By 

using the C and D solved for in the vector potential approach 
nm nm 

the equations for the S parameters and the <j> = 0, 45° and 90° co- 

and cross-polarization patterns will be derived in this section. The 

C and D are chosen since the vector potential approach is the one 
nm nm r 

used in this report. 

First we consider 2 guides each driven by the TE^q mode. Let 
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and be the incident and reflected voltages in guide 1 and 

and V 2 in guide 2. These are related as follows: 


V i = S ll V t + S 12 V 2 


V 2 = S 21 V I + S 22 V 2 


Where the S are the scattering matrix elements. 

The equivalent voltages are proportional to the amplitude of the 

incident and reflected 10 vector potential modes- so V^ + a^q^ and 
+ 2 

V 2 * A 10 ' re f^- ecte ^ vector potential mode is proportional to 

A^ 1 plus that contributed by the aperture currents using the Greens 

dyadic derived in Sec. 2.3. The latter contributes — — C, J- . 

IP 10 

1 e 0 1 10 

Thus V a - (A + — — C in ) where the negative sign is due to 

1 iU - ,B 10 1U 

electric field's dependence on the derivative with respect to 2 of 
the A^®. Knowing this, then: 

- (A ^ + — — C 1) = S A^+S A 2 

kA 10 jB °10 ' b ll A 10 *12 10 (2.83) 

- (A ^ + — — — C 2) — s A ^ 2 

U 10 jS 10 L 10 ; b 21 A 10 + S 22 A 10 


If only the 1st guide is driven the following holds: 


21 


i c io £ 0 

8 10 A 10 


and S„ = -1 

WlO 


(2.84) 


Using this the coupled power from 1 to 2 is 20 log( I s 2 i I ) an< ^ the 
reflection coefficient in guide #1 is S^. 
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In practice the only scattering matrix parameters that are useful 
are those that relate the amplitudes of the dominant propagating modes 
in each waveguide. Therefore the scattering matrix parameters for the 
other modes that may exist in the aperture aren't included. For the 
cases where there are more than two apertures the scattering parameters 
(S^^, S 4i’ ‘ • • ) are found by exciting one guide and computing 

the amplitude of the propagating modes excited in the other guides 
and then solving for the S one at a time. The procedure is a 
straightforward extension of that described above for obtaining 

The next topic in this discussion is the derivation of the far 
zone radiation patterns of the array. We use the fact that the far 
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/ ,, nm i 
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y )e x y dy dx 


( 2 . 86 ) 


J 

m,v 


rr 


V'd 1 $ (x 1 ,v’ ) e 
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(jk x' + jk y’ ) . , , 

J x J y ay dx 


0 J 0 n > m 


(2. 87) 


where k 

x 


k sin 9 cos <p and k = k. sin 0 sin $> 
o y 0 


Once chese two have been solved they can be substituted into the 
following two equations to find E^ and Eg in the far zone region. 
(Stutzman and Thiele, 1981, Sec. 8.1): 


E 

— <f> 


4—^ e^O r (J cos <J> - J sin 4> ) a 
2 irr m,y r m x —0 


_jk -jk 0 r 
2 rr 6 


( 2 . 88 ) 


(J cos d> + J sin<J>) a cos0 
m,x m,y —9 


(2.89) 


F or the multi aperture problem the J m „ and y terms need to be changed 
to account for the added current terms. 

Terms due to the other aperture currents, shifted in phase due to 
their physical spacing must be included. This gives for (2.86) and 


(2.37) : 


fj 


j (k x’ + k y 1 ) 
e x y' 


[Ci ip (x\ y’ ) + 
nm nm 


C ‘ * (x\ y» ) e 

vim nm 


0 •'o 

j ( k D + k D ) 
xx y y +> 


. ] dy' dx 


(2.90) 
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m.y 



j (k x' +k y' ) 

e X y [ D 1 
nra 


$ (x',y' ) + 

nm 



nm 


j (k D + k D ) 

0 (x\v’) e x x y y +. . . ] dy>d x ' 

nm 


(2.91) 


where the spacings are all referenced to guide #1. These, when added 
to (2.85) and (2.86) give the far zone field patterns for the given 
array . 

The co-and cross-polarization fields as defined by Ludwig are the 
fields measured by a horn polarized in the same direction as the rad- 
iating antenna for the former and as the fields measured by a horn 
polarized at right angles to the source for the latter (Ludwig, 1973). 
The co- and cross-polarized patterns provide an important set of pat- 
terns for use in analyzing the cross polarization interference between 
adjacent scanning areas. The definitions of the two are: 

R(G,<t>) = E ( 0 , ) • (sin <p a + cos pa) (2.92) 

— — 6 — <p 

C(9,<t>) = E(0,p) • (cos <t> a „ - sindi a ) 

— — 9 — <p 


where R(0,P) is the reference or co-polarization field and C(g,p) the 
cross-polarization pattern (Ludwig, 1973). Using this the resulting 
pattern equations are: 
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-j 

m.x 


=0 


J 

m,y 


cos e 


(2.93) 
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R(6,$)| = 45 = = (1-COS0 )- (l+cos 8) (2.94) 


cce,*) | ^ = — y -~ ( l+cos 0) + — ^’ - X (cos 9 -1) 


R(84) L _ qn o = -J 

<P - 90 m,x 


C( 9, (*>) | <() __ 0 J cos 9 

= 90 = m,y 


where Che common 


ik Q e _Jk 0 r 


factor was not included in the equations. 


With the equations now defined, the only problem remaining is the 

solution for the J and J terms. 

m,x m,y 

Using the following: 

/ . . ax, e ax (asinbx - b cos bx) 

sin(bx) e dx = 2 2 

a + b 


cos (bx) e ax dx 


e ‘ (a cos bx + b sin bx) 

2 ^ K 2 

a + b 


(2.86) and (2.87) can be solved to yield: 


V'' r ( £ ° ne ° m N ^2 f rnr [1-cos (nn) e 2k x 3 ] jk 
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(2.96) 


J = ^ ' n f WL ‘ u ‘“ ') 

m,y ' 'nm ' ab 

n,m 


on'om, 2 [ mir [1-cos (mm) e J v 




b [ (dajr) 2 


[e 2k ^ a cos (nr) — 1 ] 1 
[(nm ) 2 - k 2 ] J 

T x 


(2. 97) 
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The only problem with the use of these equations ((2.9b) and (2.97)) 

is the case where k = k = 0 (9 = 0°). At that point J =0 

y x r m,x 

except when m = 0 also . then J = tt. So this situation 

1 m,x 0 

must be tested for and when it arises, the result, found by using 
l'Hospital's rule substituted in it's place. 

That concludes the pattern analysis. Since the 2 sets of terms 
are the same for each current (different modes in different guides) 
except for the ,Dy phase terms the patterns can easily be calculated 
as a sum of these quantities. 

These two sets of parameters then conclude the analysis of an 
antenna array with mutual coupling. The S parameters yield coupled 
power and reflection coefficients while the patterns calculations yield 
the side lobe levels, cross-polarization and scan affects. 

2. 7 First Order Coupling 

Since the computer time needed to solve large arrays can be con- 
siderable, future work in ways to approximate the problem is needed. 

One possible technique is the idea of first order coupling. 

First order coupling solves for the coupling from a driven guide 
to it's surrounding undriven guides and ignores any fields that might 
be coupled back due to the undriven guides re-radiating the fields. 

The possible use of this technique could reduce computational costs 
significantly, but needs further study before it can be considered 
complete. 

The main idea behind this technique is that the magnitude of 
a field that is coupled to a neighboring guide from a driven guide 
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and then re-radiated back again should be small enough to be ignored 
Figure 2.7 may help explain it better. For the two aperture case 
(guides 1 and 2) only the is considered. The Y-^ is ignored 



Figure 2.7 Coupling in a Three Aperture Array. 


For the case where only one aperture is driven the results 
taken using the computer model for 2 and 3 apertures seem to support 
this theory (see Chapter 3). It was also found that for the 3 aper- 
ture case with the first one driven, that the coupling from the 
second to the third aperture must also be included to accurately 
predict the field in the third. Looking at Figure 2.7 this means 
that the Y , Y and Y^ must be accounted for. The reason for 
including Y is that (3) is closer to (2) than (1). So the re- 
radiated field from (2) is large enough compared to the field from 
(1) that it cannot be ignored in the case presented here where all 
three apertures are in a straight line. 

The incorporation of this idea into the computer model is very 
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easy. Basically all the terms for coupling from (2) to (1), (3) to 
(1) and (3) to (2) are set to zero. 

This is shown below: 


~ Y 11 

Y 12 

Y 13 ” 


‘ Y 11 

0 

1 

o 

Y 12 

Y 22 

Y 23 


Y 12 

Y 22 

o ! 

Y 13 

Y 23 

Y 33 


_ Y !3 

Y 23 

Y 33 _ 


Full Coupling First Order Coupling 


(reciprocity is used in these matrices). This is for the case of a 
single driven element. 

The advantage of this approach is that for large arrays with 
only iJ^q modes in the current approximation, the resulting matrix 
is a lower triangular matrix, like the one shown above. These types 
of matrices are very easy to solve using back substitution (Cheney 
and Kincaid, 1980), which does not require a matrix inversion. Also, 
the number of needed storage locations is cut almost in half. So 
a given computer could handle larger arrays. Using these results 
then, the S parameters and the radiation pattern of a single 
radiating element in a large array could be found quite rapidly, 
even if more than a single mode is used the computation time is 
still reduced in large arrays because the number of admittance 
terms to solve for is smaller. 

To use this idea in a large array with several guides driven 
is not as easy to do. This is a part that needs further study. 
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A possible technique may be to solve for the currents in a few rep- 
resentative cases (guides on the edge, on a corner and in the middle) 
using this approximation in smaller arrays. Then applying super- 
position, the entire array solution could be found. A second and 
maybe more promising idea would start by breaking the large array 
up into a smaller one. The coupling terms from this smaller array 
could then be substituted into the complete larger matrix which would 
then be inverted. This says that if first order coupling is all that's 
needed then coupling from far away guides (several rows away) and 
from re-radiated fields can be ignored. This idea would consider- 
ably reduce the time needed to calculate the admittance matrix. 

Further tests on this subject need to be run before it can be 
considered a valid approach. Some tests for 2 and 3 aperture arrays 
are presented in the next chapter. 

2 . 8 Conclusion 

In this chapter four techniques for solving for the mutual 
coupling between open ended waveguide elements in an array were pre- 
sented. They were (1) Fourier transform method, which solves directly 
for field amplitudes in the transform domain, (2) the TE and TM mode 
and (3) sinusoidal-piecewise aperture field expansion methods which 
use the method of moments to solve for the amplitude of equivalent 
magnetic currents, and (4) the vector potential approach which is similar 
to (2) except in the choice of basis functions and the manner in which 
the integral equations are derived. 
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In looking at the various approaches many similarities between 
them appear. The Fourier transform method used by Bailey (1972) is 
very similar to the TE and TM modal function approach. 

They both solve for the admittance matrix that relates the inputs 
(the source modes of the driven guides) to the resulting fields which 
are approximated by TE and TM modal functions. The latter uses the 
method of moments to derive the admittance equations and the former 
solves 1 the reaction integral for admittance directly* using Fourier 
trans forms . 

The three method of moments approaches set up the same problem, 
but differ in choice of basis functions. Harrington and Luzwick and 
the vector potential approach use modal functions to reduce the size 
of the resulting matrix. The sinusoidal-piecewise basis functions 
used by Fenn et al, on the other hand, are used because the resulting 
integral equations are easier to evaluate. Another advantage of this 
approach is that for large numbers of basis functions per aperture 
the guide edge singularities can be analyzed (see Fenn et al, 1982). 
The main drawback with it though is that for square guides large 
numbers of patches are needed (24 was found to be necessary) per 
aperture. This results in very large matrices that must be solved 
for and inverted, preventing the solution of large arrays. 

The vector potential approach was the one chosen for use in 
the computer model presented in this report and was derived in 
greater detail. The reason for it's choice was that by using the 
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vector potential modal functions the integral equations were more 
straight-forward to work with and solve. It is also easier to in- 
clude several, including orthogonal modes, in the solution, and 
finally, the resulting equivalent magnetic corrents are easy to work 

with in solving for the radiation patterns. 

The rest of the chapter (the last three sections) present math- 
ematical techniques and parameter derivations that were used in 
designing the computer model along with a possible method for ap- 
proximating large arrays. 


CHAPTER 3 


COUPLING EFFECTS 

3.1 Introduction 

A computer program was written to solve the equations derived 
in the last chapter using Simpson's Rule (see Appendices A and B for 
listing and user instructions). It is a very flexible program that 
can accept any reasonable number of modes and array elements. It can 
also be configured to output first order coupling results. In this 
chapter samples of data run on this program will be presented. 

Some of it will be compared against published data to test the valid- 
ity of the model, while the rest will be used to provide more insight 
into the mutual coupling problem. 

3.2 Single Aperture 

The analysis of a single aperture in a ground plane provides 
much needed information for use in computer runs involving larger 
arrays. The results provide a means for determining the best mag- 
netic current approximation to use (which modes). The data also 
provides insight into the problem of a waveguide to free space inter- 
face by showing the types and magnitude of the fields excited there 
and the size of the reflected wave. 

The determination of the number of modes needed for a good 
approximation is the first topic to be considered. The symmetry 
of the incident TE^q E field is even about the midplane of the y axis 
and even about the midplane of the x axis. The only fields excitpd in 
the aperture then should be ones of like symmetry (Borgiotti, Nov. 


59 


60 


1968) . It follows then that the testing and basis functions should 

be ip with n = 1, 3, 5, . . . and m = 0, 2, 4, . . . and v with 
nm nm 

n = 1, 3, 5, . . . and m = 2, 4, 6, . . Table I shows the results 
of computer runs made using various combinations of these modes. The 
aperture is square and 0.6A in width. Shown are the normalized aper- 
ture E field modal amplitudes, (see Appendix C for normalization used). 
Table II shows the results for a like set of modes for a rectangular 
aperture where a = 0.8A and b = 0. 4A using the same normalization. 

It can be seen from these two tables that ip ^^ , an< i ^]_2 

appears to be the best choice. It has the fewest number of modes 
that shows no change in when additional modes are included. Also 
any added modes are of magnitudes that are 3% or less of C^q. 

It should be noted when looking at these tables that there appears 
to be a strong interaction between orthogonal modes of identical n 
and m. In Table I C ^ changed by almost 50% when was added. This 
happened for all such cases listed in both tables. This interdepen- 
dence is expected due to the fact that for a given nm waveguide mode 
where m ^ 0 there is both an x and y component to the E field. It 
is because of this phenomenon that should be included in rectan- 

gular waveguide solutions (Table II). The magnitude of is less 
than 2% of C^q but it's addition means a large change in So 

when using a ip with m 4- 0 it's corresponding $ should be included 
6 nm v 6 nm 

to insure proper results. 

After the number of modes was determined a comparison of reflec- 
tion coefficients, found using the vector potential approach, against 
ones, found by Fenn et al. using the piecewise sinusoidal functions. 
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was run. The first set of points is for a rectangular aperture vary- 
ing the width (measured in number of X ) but not the width to length 
ratio. Plots of the results are shown in Fig. 3.1. The comparison 
data is read off a graph included in the article (Fenn et al, 1982) 
for 7 sinusoidal patches all x directed. The next case was a square 
aperture, again varying the width. It is shewn in Fig. 3.2, Fenn 
et al's data were also read off a plot included in their article. 

These two sets of plots show the excellent agreement that exists 
between the two models. 

Finally the plots of the radiation patterns for a single square 
aperture are presented in Figures 3.3, 3.4, and 3.5, Each plot is 
for a given phi angle (0°, 45° or 90°) with the theta varying from 
-90 to 90°. The coordinate system is shown in Fig. 2.1. The phi = 0° 
is equivalent to an H-plane pattern and phi = 90° to an E-plane 
pattern. The phi = 45° plot shows both the co- and cross-polarized 
patterns. For <f> =90° and 0° there is no cross-pol. pattern due to the 
symmetry of the single aperture antenna. These plots are included here 
for comparison against results presented in l?tP r serf'on^. 

3.3 Two Aperture Analysis 

In this section further comparisons of computer derived results 
against published data will be presented. These will be followed by 
a look as some tests of the first order coupling hypothesis for two 
aperture arrays. 

The first comparison presented is one between results from this 
model and some of the results presented by Luzwick and Harrington 
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HS 

(Luzwick and Harrington, 1982). They presented plots of the Y 
admittance term (see Section 2.4) for coupling between two apertures 
at a constant distance but at varying angles referenced to each other. 
The comparisons are shown in Fig. 3.6, where dimensions are given in 


terms of X. The computer data used a single mode approximation 
(<J> ), similar to what was used by Harrington and Luzwick. The agree- 

ment here is very good. 

The next comparison is against data published by Mailloux 
(Mailloux, Nov. 1969). In his article, results are presented for 
coupling between two apertures that each support two orthogonal modes 
(TE^q and TE^) . The S parameters for coupling from the dominant 
TE-^q mode in the driven guide to the TE-^q and TE^ modes in the un- 
driven neighboring guide are plotted. The results are for two aper- 
tures at constant distance, but with varying angles again. Compari- 


sons of SI , 2 

TE 10 TE io 


and $ rpp 2 

TE 10 TE 01 


are shown in Fig's 3.7, 


3.8, 3.9, and 3.10, for two different spacings. Dimensions are in 


terms of X . Here again good agreement can be seen. and 

modes were used in the computer approximation. 

It should be noted that the agreement with Harrington and 
Luzwick' s . results appears to be better than the agreement with 
Mailloux's results. Harrington's plots were larger and easier to 
read data off of and that is the reason for the better fit on those 


plots . 

Now that the validity of the computer model has been established 


a few tests of the first order coupling hypothesis, whose validity 
needs to be checked, will be presented. As a reminder, the idea 
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behind this approximation technique is to assume that the fields re- 
radiated by neighboring undriven elements are small enough to be left 
out of the calculation, fora more detailed discussion see Sec. 2.8. 

To test the first order coupling idea, computer runs were made 
for two co-linear apertures mounted horizontally using the four mode 
approximation determined in the last section. Figures 3.11, 3.12, 

3.13 and 3.14 show comparisons of the mode amplitudes found using 
first order coupling against the mode amplitudes found using the 
complete matrix. It can be seen that the effects seem to be larger 
on the driven element than on the neighboring undriven one. This is 
because the effect on the undriven element due to the removal of the 
re-radiated fields is a second order effect so, it’s change will be 
small. Since it is so small the next set of plots will look only at 
the driven aperture amplitudes. 

Figures 3.15, 3.16, 3.17 and 3.18 show the comparison of first 
order mode amplitudes verses the complete matrix amplitudes for two 
different aperture arrangements. Figures 3.19 and 3.20 show the same 
set of amplitudes but for a pair of rectangular apertures. 

Finally to test the effect on the patterns, Fig. 3.21, 3.22, 
and 3.23 show a comparison of the patterns for the three configu- 
rations with square guides and 0.8 A spacings. All are for only one 
element driven. The dashed curves which are the first order coupled 
plots are almost indistinguishable from the solid curves which are the 
complete plots. Therefore these patterns can be considered virtually 


identical. 
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This concludes the two aperture analysis which helped to provide 
some validity to the computer model derived using the vector potential 
approach. Further testing of first order coupling will be presented 
in the next section on three aperture arrays, and any comment on it 
will be with-held till that time. 

3. 4 Three Aperture Arrays 

The results shown in this section will help to further test the 
first order coupling hypothesis. Computer data for three co-linear, 
horizontally mounted apertures were collected. These data show that 
the first order coupling approximation does give valid results. Also 
included in this section will be a few plots of the patterns of three 
different arrays where the beam is being scanned. These patterns 
illustrate the beam forming capability of the computer model. 

The first order coupling tests performed were similar to the 
ones presented in the previous section but with three apertures. In 
this case the co-linear apertures were mounted horizontally at various 
spacings. The first set of plots is for the left most anerture 
driven. The 4 mode approximation discussed in Sec. 3.1 ^12 ’ 

'^q and $-^ 2 ^ is used here for all the data in this section. Figures 
3.24-3.29 show the mode amplitude results for all three apertures of 
this test. Figures 3.24 and 3.25 show the mode amplitudes of the 
driven guide while Fig's. 3.26-3.29 show the same for the middle and 
right most guides. In this situation only, coupling between the 
second and third apertures was included. Preliminary tests showed 
that without it large errors in the third aperture's field resulted. 


c-2- 
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The reason for this is that the coupling from the driven to the third 
is small enough that the re-radiated fields from the second aperture 
could have a sizable affect on it. 

The next set of data is for the same aperture arrangement, but 
with the middle guide driven. Shown in Fig's. 3.30-3.31 are the mode 
amplitudes of the driyen guide. With no coupling between the end ele- 
ments their amplitudes should be similar to the two aperture case 
of Fig's. 3.13 and 3.14 and will not be repeated here. 

Finally, Fig's. 3.32 and 3.33 shows the comparisons of the radia- 
tion patterns for the two cases described above. Each shows the full 
coupling patterns and the first order coupling patterns for three 
aperture arrays with 0.7A spacings. There appears to be little 
difference hetween the first order and the full coupling plots. This 
is an expected result due to the near identical patterns found in the 
two aperture case, Fig's. 3.22 and 3.23. 

This concludes the computer tests for the first order coupling 
idea. Some experimental tests were performed and are presented in 
the next chapter. These experimental results agree with the computer 
data presented so far. The preliminary conclusion to be made from 
this testing is that the theory gives a fairly accurate solution to 
the mutual coupling problem, at least in the cases presented here. 

The largest difference found was less than 5% for the magnitudes and 
a few degrees for the phase. It was also found to have very lictle 
affect on the dominant TE^q mode, which is the largest mode present. 
These encouraging results warrent further testing of the theory and 
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more work on formulating a viable way to apply it to the types of 
antenna problems encountered today. 

The final set of plots in this section are examples of the beam 
scanning properties of a three element array. The proper phase shifts 
in the driving modes of each aperture were entered to give 10° and 
20° scans for a colinear array of three apertures at three different 
spacings. The results are shown in Fig's. 3.34 to 3.39 (again the 
4-mode solution was used here). The plots are for <j> = 45° and ® = 0° 
(each normalized by it's own largest value). Since there is only a 
single row of apertures the change in the shape of the patterns for the 
<f> = 90° was slight. There was an overall reduction in level, but the 
shape didn't change with change in scan angle. So these patterns will 
not be presented here. 

Before leaving scanned beams a look at the patterns show some 
interesting, even though expected results. First it's easy to see the 
super- imposed single aperture radiation patterns. These are products 
of the array factor, and can be seen in Fig's. 3.4, 3.5, and 3.6. The 
characteristic three aperture array factor can also be seen with it's 
single side lobe between main and grating lobes. Due to the large spac 
ings necessary in an aperture antenna the grating lobe becomes sizable 
with very small scans of the beam. Therefore the resulting patterns 
do exhibit the general behavior of array antennas. 

This concludes the three aperture results and the first order 
coupling computer testing. As was stated, the first order theorv 
appears to provide an adequate approximation to the mutual coupling 


and should be studied further. 
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3.5 Modest Sized Arrays 

The results to be presented in this section are patterns of 
modestly sized arrays. The purpose is to compare the ideal array 
patterns with ones that include mutual coupling. These comparisons 
will help to show the necessity of including mutual coupling in larger 
arrays. A final note on edge effects will complete the discussion. 

The first set of plots is the radation patterns of a 3 x 3 array 
of square apertures (0.6A x 0.6A) in a rectangular lattice with 0.7A 
row and column spacings. The plots show a comparison of the ideally 
derived patterns (no coupling present) with the patterns found using 
the mutual coupling. These are shown in Fig's. 3.40, 3.41 and 3.42. 
Figures 3.43, 3.44 and 3.45 are the patterns for a 5 x 5 array of 
square apertures with like spacings. The 7x7 array case is shown 
in Fig's. 3.46, 3.47 and 3.48. For the 3x3 and 5x5 array cases, 
a 5 mode solution (^^g> ^ 12 ' ^ 30 ’ ^01 anc ^ WaS usec *' E° r t ^ ie 

7x7 array only, a single mode, the dominant TE^q mode was used 
due to the limited matrix size of the program. All of the patterns 
are for the case of all the apertures equally driven with equal phases 
and amplitudes. 

In looking at these plots the first conclusion to be drawn is 
that the coupling between columns must be smaller than the coupling 
between rows. This is seen in the fact that all of the 4 1 = 0° (H- 
plane) comparisons show little difference between the patterns wicn 
and without coupling. There is a slight difference, but it didn't 
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always show up when the two sets of plots were superimposed. In look- 
ing at the data presented in Sec. 3.3 at Figures such as 3.7 or 3.9 
it can in fact, be seen that coupling is smaller into an aperture next 
to the driven element then into one directly above it. So these 
patterns bear that out. 

It can also be seen in these figures that the effect on the co- 
polarized main lobe and the first side lobe (in the larger cases) seem 
to be small. It wasn't until the second and higher side lobes that 
there seemed to be a difference. So the mutual coupling affects seem 
to be greater on the sidelobes that are at modest angles from the main 
lobe. The effect on the cross-polarized signal (p = 45°) though was 
considerable. The coupling not only caused amplitude changes, but 
also shifted the position of the nulls. This result is due to the 
fact that the mutual coupling terms allow for additional energy in 
orthogonally directed fields to be present in the apertures that are 
cross-polarized in nature (of different symmetry than the incident 
fields). This additional energy, unaccounted for in the ideal case, 
then caused the large changes. It also seems that the differences in 
the 7x7 case are small. This is due to the fact that only one mode 
was used in the computer solution , thus not accurately predicting all 
of the coupling present. This shows the need for including more than 
just the dominant TE^q mode in the aperture field approximations. 

The next comparison is that of the E-plane (p = 90° ) pattern of 
a single aperture radiating in various sized arrays, where all other 
elements are terminated in matching loads but not driven. Figure 3.^9 
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shows these patterns for three array sizes. The first is that of 
a single element array. The other two are for a 5 x 5 and an 11 x 11 
array where the center element is the driven aperture. These are all 
square apertures with the arrays in a rectangular lattice with 0.7A 
row and column spacings. For the single and 5x5 array case, the 5 
mode solution discussed above was used while a single mode was used 
in the 11 x 11 array case. 

f 

It is easy to see the effects that mutual coupling has on the 
patterns as more apertures are added. These additional apertures 
absorb energy that would normally be present in the pattern and 
scatter energy that then interacts with the original radiated fields. 
These two effects then account for the changes seen. It is interest- 
ing to note that as more elements are added, the pattern gets slightly 
more directive. 

The analysis of the E-plane patterns of the center element of 
a large array was also looked into by Amitav, Galindo and Wu (Amitay, 
Galindo and Wu, 1972). They experimentally measured this pattern for 
the center element of a 19 x 19 array of square apertures. The array 
had slightly closer spacings than that used above. Their results 
show that at approximately 30° and 70° a sharp null appears in the 
pattern. The beginnings of this null can be seen in the 11 x 11 plot 
shown in the figure. While not as sharp, there are dips in the 
pattern chat are 4 dB oelow cne single aperture pattern. So it seems 
that the computer derived patterns do in fact give a fairly accurate 


picture of tne true patterns. 
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Turning to the topic of edge effects, it is necessary to look 
at the early work on mutual coupling. In the early analysis the 
finite arrays were approximated by arrays with an infinite number of 
rows and columns. This allowed the radiated fields to be represented 
by Floquet modes and was thus easier to analyze. (see Amitay, 

Galindo and Wu, 1972 for complete analysis). The problem with this 
type of solution is that it predicts the effects on central elements 
in a large array fairly well, but it doesn't give accurate results for 
elements closer to the edge. This difference between coupling affects 
of center verses edge elements is known as the edge effect. It can 
be explained intuitively by considering the neighboring apertures that 
an element "sees”. Center elements are surrounded on all sides by 
apertures while edge elements "see" fewer closely spaced aper- 
tures. To illustrate this point the patterns of the center element 
and the upper right hand corner element of a 5 x 5 arrays are compared 
in Fig's. 3.50 and 3.51. The elements are square (0.6A x 0.6A) with 
.7A row and column spacings in a rectangular lattice. A single mode 
approximation was used here. From these plots it can be seen that 
the edge element's patterns are shifted to one side due to the fact 
that it doesn't have neighboring elements on two of it's sides. 

So it can be concluded that due to the edge effects, the finite 
array type of analysis performed here is necessary to insure proper 


results . 
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Comparison of Patterns of a 3 x 3 Array; ip = 0°. 
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Figure 3.42 Comparison of Patterns of a 3x3 Array, 4» = 90 
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Figure 3.4& Comparison of Patterns of a 5 x 3 Array, <j> = 45 
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Figure 3.45 Comparison 
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Figure 3.46 


Comparison of Patterns of a 7 x 7 Array, 4> = 0 
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Figure 3.47 Comparison of Patterns of a 7 x 7 Array, <j> - 45 
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1 - Single Element Array 
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3 - 11 x 11 Array 
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Figure 3.49 


Comparison of E-plane Co-polarized Patterns for Center 
Element of Three Different Sized Arravs. 
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3.6 Conclusion 

In this chapter various computer derived data, found using the 
vector potential approach, were presented. These results helped to 
shed some light on the mutual coupling problem by looking at it 
from several points of view. 

First the single aperture array data was used to determine the 
proper approximation to use in the later analysis. After this, the 
validity of the computer model was established by comparing the 
results with those of several published articles. These comparisons 
along with the data to be presented in the next chapter show that 
the model is in fact reasonably accurate. 

Some preliminary testing of the first order coupling idea was 
then presented. The two and three aperture tests showed that the 
first order coupling approximation did yield valid results. The 
testing, while incomplete, was sufficient to warrant further study 
on this idea. 

The rest of the chapter dealt with various sized arrays and how 
mutual coupling affected their performance. This included a look 
at the beam scanning properties of small arrays and concluded with 
a discussion on edge effects. 

Finally, a word regarding the radiation patterns presented here. 

The patterns were each normalized by their own largest value. That means 
that the <j> = 90°, <p = 0° and <j>= 45° patterns were each normalized by their 
own peak values- The co- and cross-polarized patterns for + = 43’ were 
both normalized by the same value. This shouldn't cause any problems 
for the non-scanned plots as the peaks for all three coincide at 9=0. 
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In the scanned plots of Fig.'s 3-34-3.39 though this is not the case. 
If the i p= 45° and <p = 0° plots of each figure were normalized by the 
same value, the overall level of the <f>= 45° plots would be slightly 
lower. The shapes of the patterns wouldn't change though, so the effect 
of scanning the beams is still properly illustrated by these plots. 


CHAPTER 4 


EXPERIMENTAL VERIFICATION 

An experiment was conducted to test the accuracy of the computer 
derived results. The tests measured the coupling C S 21 ) between two 
apertures in various configurations on a large copper plate. The 
measurements were made with a Hewlett Packard 8409B Automatic Network 
Analyzer with version 118b 3D software. The test stand was a 2 ft. 
square ground plate with the waveguide fed apertures radiating to a 
room covered with absorber material. 

To facilitate ease of manufacture and to reduce set-up times, a 
general purpose test stand was constructed. A 1/4" copper plate 
2 feet square was milled so as to allow 5 X-band waveguides to be 
soldered directly to it. The hole configuration is shown in Fig. 4.1. 
The spacings were chosen by taking into consideration the waveguide 
flange sizes and the amount of material necessary to assure a strong 
plate that wouldn't warp during soldering. 


T 

.034 
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Figure 4.1 Test Set-up. Distance in Meters. 
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Two tests were run on this set-up. The first measured the 
coupling between two apertures at various spacings and angles in terms 
of for direct comparison against the computer data. The second 

was a test of first order coupling. The and was measured 

between two apertures with a third aperture present. This third 
aperture's waveguide was terminated in a matched load. These measured 
scattering parameters were then compared against the parameters 
measured for the coupling with no third aperture present. The differ- 
ence then would be the re-radiated third aperture fields and thus the 
relative effects of the re-radiated fields could be determined. 

In all these tests the unused apertures in the plate were covered 
with foil, so as not to interfere with the coupling between the two 
apertures being tested. For example to test the coupling between two 
apertures side by side at .034m spacing, the analyzer would be con- 
nected to the guides for apertures 1 and 2. Apertures 3, 4, and 
5 would then be covered with foil. 

The first set of data is shown in Fig's. 4.2 and 4.3. It com- 
pares the experimental results for two co-linear apertures, mounted 
side by side at 3 different spacings, against data from the computer. 
The computer results for the tests in this section used the four mode 
approximation discussed in Sec. 3.2 ^ 12 ’ '^30 anc * ^12^’ Figures 

4.4, 4.5, 4.6 and 4.7 show the same comparison but for two other 
apercure configurations. In these tests the frequency was scanned 
from 8 to 12 GHz to give a set of data points for each arrangement. 
By aoing this, the size of the apertures and their spacings would in- 
crease in terms of wavelength with each increase in frequency. 
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The test data seemed to deviate from the computer results most 
at 9-10 GHz. At first it was thought that resonances in the set-up 
might be responsible, but after a more complete run was made, the true 
cause became apparent. Figures 4.8 and 4.9 show the $ 2 ^ and S ^ mag- 
nitude plots taken directly off the network analyzer for runs with 
very small spacing between data points. The experimental curves are 
seen to oscillate about the drawn in computer derived lines. The ver- 
tical axis is dB of loss so they are in fact negative quantities. 

The cause of the oscillations is probably due to the irregular- 
ities in the system. One possible reason being the fact that the 
ground plane wasn't infinite. This could account for part of the 
oscillations since the amount of interference would depend on the 
distance between the source and the diffracting edge in wavelengths. 
This would then vary as the testing frequency is scanned. A second 
source of error is the test set-up itself. It was calibrated for a 
given location of the cables and their connections. These cables 
are then moved and connected to the test stand, destroying the calib- 
ration. This would explain the fact that the plots in Fig's. 4.8 and 
4.9 give slightly different numbers then the data shown in Fig's. 
4.2-4. 7 at specific frequencies. While the average line for the 
experimental data was always about the same the overriding noise 
seemed to move up and down this line, depending on the day the 
test was run. Tests run on the same day with the same calibration 
didn't show this variation. If the noise could be averaged or filtered 
out of the experimental data, it would produce a set of curves that 
appears to be close to the computer derived data. So it seems possible 
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at this point to say that the tests agree with the computer data within 
the limits of the system noise. 

The second test was performed for two co-linear apertures with 
a third undriven one present. The result for two different aperture 
arrangements are shown in Fig's. 4.10 and 4.11. The data measured 
with a third aperture present is compared against that found in the 
previous test. These plots are consistent with the first order 
computer tests presented in the last chapter. The presence of the 
third aperture seemed to have little affect on the and S in 

Fig. 4.10, while the variation in S in Fig. 4.11 could be explained 
by the irregularities seen at 9-10 GHz discussed above. 

The final measurements made tested the level of the analyzer's 
internal noise. It was measured by terminating the analyzer's out- 
puts in matched loads and measuring the amount of signal present. 

The tests showed the worst case noise to be at approximately -60dB. 

This would explain the large differences seen in the d = .063 side- 
by-side test (Fig's. 4.2 and 4.3). It's magnitudes were approaching 
the point where the analyzer noise was large enough to have an effect. 
Also tested was the amount of reflected signal from the absorber in 
the measuring room. It showed that the reflection was buried in the 
level of the analyzer noise. So these two noise sources are small 
enough not to have much effect on the results. 

This concludes the experimental analysis of mutual coupling be- 
tween apertures in a ground plane. The results for the most oart are 
consistent with the data found by the computer using a 4 mode solution. 
While the test results were slightly irregular, they were sufficient 
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enough to verify the accuracy of the analytical model. 
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CHAPTER 5 


CONCLUSIONS AND RECOMMENDATIONS 

5.1 Conclusion 

The work described dealt with the analysis of mutual coupling in 
finite arrays of rectangular apertures. The analysis started with a 
discussion of the theoretical nature of the problem. It followed that 
with a look at the coupling effects in various arrays, which included 
comparisons against both published and experimental results. The con- 
clusion to be drawn from this work is that the analysis and approxima- 
tions used do in fact vield accurate results for the mutual coupling 
problem. 

The various methods used to analyze the mutual coupling in arrays 
were discussed. From this discussion it was concluded that the vector 
potential approach was the best one to use for the type of array being 
analyzed. The required mathematical manipulations used in deriving 
the computer program are also included. 

Using the computer programs, results were obtained for various 
sized arrays. The single element array data determined the best modal 
approximations to use. The two aperture case compared the computer 
data against published data. The other configurations looked at pre- 
sented results that showed the effects of mutual coupling on larger 
arrays. It also looked at beam scanning in the three aperture antennas. 
Included with all this analysis was a considerable amount cf data for 
testing the first order coupling theory. This data showed that the 
first order coupling theory could be used to approximate the coupling 
in arrays fairly accurately. It was included here for use in future 
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studies on antenna arrays in the hopes of reducing the computational 
time needed to solve large arrays. 

Finally, in the last chapter the results from an experiment that 
was carried out was compared against computer derived data. While the 
experimental data was affected by system irregularities it did seem to 
show that the computer-derived results were correct. 

5.2 Recommendations 

Future study on the mutual coupling problem should: 

1) include the non-matched source to waveguide case. 

2) study arrays of circular apertures. 

3) look at flared horns, both circular and rectangular. 

4) consider the no-ground plane case. 

5) look at possibly formulating a viable technique that uses 
the first order coupling. 

The non-matched load case is one that can be easily included in 
the current computer model. If the waveguide is l meters long and the 
guide to source interface has a known reflection coefficient (T) , then 
to include it's affect a factor (1 + TTe ^ 1024^ must b e multiplied 
by the Green's Function term for the TE^q P art the waveguide field 

(G xx for nm = 10 in eqn. 2.36). In this factor, T is the total am- 
plitude of the dominant mode excited in the guide due to reflection 
at the aperture and coupling from all the other guides. Once this 
has been done, tests could be run to see if the effect of this on the 
system results is large enough to warrant the inclusion of this type 
of mismatch in future studies. 

The other problems are much more involved and would require con- 
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siderable amounts of work. The no-ground plane analysis should be 
conducted to determine if the affects of removing it are large enough 
to require that they be accounted for in antennas where the ground 
plane isn't present. If the effects are small enough, then the easier 
to solve infinite ground plane case could be used as an acceptable 
model. The circular aperture and flared horn elements in arrays should 
also be analyzed since these are also common array elements in use 
today. The final topic, the first order coupling theory, is one that 
needs further computer testing and much theoretical analysis before 
it can be considered complete. This analysis should be performed 
since it may lead to more efficient numerical procedures for analyzing 
mutual coupling effects in large arrays. 
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APPENDIX A 

PROGRAM LISTING 


original 

* POOR 


PAGE is 
QUALITY 


10 
20 
3 C 
40 
50 
60 
70 
80 
90 
100 
no 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
2b0 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 
460 
470 
480 
490 
500 
510 
520 
530 
540 
550 
560 
570 
580 
590 
600 
610 
b 2 0 
630 
640 
650 
660 
670 
b 8 0 
690 
700 
710 
720 
730 
740 
750 
760 
770 
780 
790 
800 
810 


C 

c 

c 


c 


1 0 
20 
30 
40 

C 


50 


60 

C 


C 


70 

C 


80 


SOLUTION OF MUTUAL COUPLING IN ANTENNA ARRAYS 
Of RECTANGULAR APERTURES. MATRIX Of 360X360 
ALLOxEO. 

REAL A, b, LA, LAM, PI ,EO,TP,E,REM,XAt 180) , IER,RP,DX ,OY,X, Ih 
REAL D,EF,S,DXP(180) ,DYP(180) ,DT,LB,DELl ,0EL2 
REAL X(91),Y2(91),Y3(91),MX1,MX2,KT(360,4),MA3,MX4 
REAL PP(lS0),PHR,iDjPR,PHG,D3,Y4(91),Y5(91) 

REAL Y1 1 (91 ) , Y 12(91 ) , T13C91) ,Y 14(91 ) ,MXT 
COMPLEX BT, BETA ,08 LINT, Z1 ,Z2, 23,24, Yl , nl ,H2 ,11,12,13,14 
COMPLEX 15, I b ,C,G,EM1 , EM 2 , J MX , J M Y , EM 1 , EM4 , R C , V 1 , NORM 1 , V 1 5 
COMPLEX* 1 6 AM(3b0,3e0),V(3e0),NV(3oU,3b0),VI(360) 

CHARACTER * 1 ANS , TEST , TEST2 

INTEGER CM(6,2J.P,U.N,M,ND,NA,NC,NT,NR,NAP,VP(180),RTP 
EXTERNAL TP,Yl,Hl,ri2,Il,I2.I3,I4,I5,lo,RP,G 
COMMON P,0,N,M,A,e,K,£,D,DX.UY 
SET UP N V AS AN IDENTITY MATRIX 
DO 10 M 1 9= 1 ,60 

NVlM19,M19)s( 1 .0,0) 

CONTINUE 
FORMAT C A 1 ) 

FORMAT (IX, 'REFLECT COEFF OF GU I OE I • , 1 X , 1 2 , 1 X , ' 1 S : • ) 

FORMAT C2X, 'INPUT THE » OF APP. IN RGx » ' , 2X , 12, XX, 

1 'AND THE KUx OFFSET' ) 

SET UP IRE NElDED CONSTANTS 

TYPE*, 'SOLVES MULT I -APPERTUR E PROBLEM • 

TYPE*,' USING THE VECTOR POTENTIAL APPROACH' 

TYPE*, 'INPUT THE » PSI AND PHI FUNCTIONS' 

ACCEPT*, NC.NL 
NT=NC+NO 
RTP = 0 

DO 50 J 22= 1 , NC 

TYPE*, 'INPUT THE ' , J22, ' PSI CGEFF.!' 

ACCEPT*, CM(J22,1),CM(J22, 2) 

CONTINUE 
DO GO J23=l , NU 

TYPE*, 'INPUT THE' ,J23, 'Phi COEFF.:' 

ACCEPT*, CM (NC*j23, 1 ) ,CM(nC + J2 3,2 ) 

CONTINUE 

TYPE*, 'INPUT A/ LAMBDA , 8/ LAMbDA AND LA MbU A ' 

ACCEPT*, LA, Lb, LAM 

CM=THE PHI oPSI COCFFS ND=*PH1 nC=»PSI NT=TUTAL» 

A=LA*LAM 

8=Lb*LAM 

TEST=' Y' 

PI=3. 1415926 
E0=8.854E-12 
AL=(A*.001)**2 

X=2*P1/LAM 

REM = 4.0*(b*LUG( ( A + SOR T ( A * * 2 + B* * 2 ) )/b)+A*LUG( (b*SUHT(A**2 
1 +B**2))/A)) 

NA = 0 

TYPE*, 'ALL SPACINGS ARE TO BE GIVEN IN FRACTIONS OF LAMbDA' 
TYPE*, 'INPUT THE • OF ROxS AND SPACING BETxEEN THE RUxS’ 
ACCEPT*, NR, DELI 

TYPE*, 'INPUT THE SPACING BETXEEN ELEMENTS XITH1N A ROx ' 

ACCEPT*, UT 

OT=DT*LAM 

NRsiROxS DT=ELEMENT SPACING DELlsRUX SPACING 
NAP£=0 

TYPE*, 'REGULAR LATTICE ARRAY? CY/N)' 

ACCEPT20.TEST2 
IF (TEST2.EQ. ' Y ' ) THEN 

TYPE*, 'INPUT THE EVEN RUXS I OF GUIDES & KUx UFFSET ' 
ACCEPT*, NAPE, DEL2E 

TYPE*, 'INPUT THE ODD RUxS » OF GUIDES & ROx OFFSET' 
ACCEPT*, NAP0,0EL20 
DO 9 U U50=l , NR 

IF ( MOD (JbO,2).EQ.U) THEN 
DO 70 J55=l , NAPE 

DXP(U55+NA)=( J55-1 )*DT+DEL2£*LAM 
DYP(u55+NA)=CJ50-1)* LAM* DELI 
CONTINUE 
NA=N A+NAPE 
ELSE 

NAPE,DEL2£=»iSPACING IN EVEN ROxS NAP0,DEL20=0DD ROxS 
DU 80 JbO=l , N APO 

OXP(UbO+NA)=( JbO-1 ) * DT+ DEL 20* L AM 
DYP(Jb0*NA)=(J50-l)*LAM»DELl 
CONTINUE 
NA=NA+NAPO 
ENO IF 
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PAGE 2 


620 

90 

CONTINUE 

8 30 


ELSE 

840 


00 110 J30=1,NR 

850 


TYPE40 , j30 

8 6 0 


ACCEPT* , NAP, DEL2 

870 

880 


00 100 J31=1.NAP 

0XPCNA + J31 )=DEL2*UAM+DT* ( J 3 1 ” 1 ) 

890 


DYP(NAfJ3l)=DELl*LAM*(J30-l) 

900 

100 

CONTINUE 

910 


N A = N A + N A P 

920 

C 

NA=TUTAL» APERTURES 

930 

110 

CONTINUE 

940 


END If 

950 


UU 120 J75=1,NA 

960 


TYPE*, 1 OX , D Y ' ,J75,DXP(J75) , D YP ( J 7 5 ) 

970 

120 

CONTINUE 

980 

C 

ThE APPERTURt ROxS 

990 


OG 310 K1=1,NA 

1000 

C 

ThE APPERTURE COLUMNS 

1010 

1020 

1030 

1040 

1050 

1060 

1070 

1080 

1090 

1100 

1110 

1120 

C 

UU 300 K2=M,»A 

If IM.EU.K2) THEN 
If CK1.EU.1) THEN 
CO 150 J=1,NC 
P = CM CJ, 1 ) 
U=CMCJ,2) 

THE COLUMNS WITH PSI0S 

UO 130 Jl=l . NC 
N=CM(J1 , 1 ) 

M = CM l J1 ,2) 

E = EF (P,U,N,M) 
b£TA=BT(N) 

1130 


0=1.0 

1140 

1150 

1160 

1170 

1180 

1190 


23=(0,”1.0)*(K**2-(N*PI/A)**2)/bETA 

21=DbLInTC0,A,0,8,Y1) 

Z2=0bLlNTCU,A,U,b,hl ) 

Z2=Z2*P*E/(2*A*a) 

If CP.EO.N .AND. U.EU.MJ THEN 

Z1 = (K**2-(P*PI/A)«*2)*U1+PEM)/(2*PI)+Z3 

1200 


ELSE 

1210 


Z1=Z1*CK**2-CP*P1/A)**2)/(2*PI) 

1220 


END If 

1230 


AM C J , J1 )=Z1*Z2 

1240 

130 

CONTINUE 

1250 

c 

THE PHI COLUMNS 

1260 


UO 140 K 4= 1 , ND 

1270 


N=CM (NC+K4, 1 ) 

1260 


M=CM(NC+K4,2) 

1290 


E=EF CP,U,N,M) 

1 300 


U=” 1 .0 

1310 


bETA=BTlN) 

1320 

1330 


Z1=U8LINTC0, A.O.b.Yl ) 
Z3=lO,1.0)*N*M*Pl**2/(8ETA*A*b) 

1340 


IF (P.EQ.N .ANU. CuEu.M) TnEN 

1350 


Zl=(-P*u*PI/C2*A»b))«(Zl+REM)+Z3 

1 3bO 


ELSE 

1370 


Z1 =-P«u*PI*Zl/ 12*A*B) 

1380 


END IF 

1 390 
1400 


Z2=dBLIN 1(U,A,0,B,H2) 
Z2=E*P*Z2/ C2*A*b) 

1410 


AM ( J , NC+K4)=Z1 t 22 

1 420 

140 

CONTINUE 

1430 

150 

CGNTINUE 

1 440 

C 

THE PHI TESTING FUNCTION ROMS 

1 450 
1460 


DO 160 K 6= 1 . ND 
P=CM CNC + K6 , 1 ) 

1470 


0=CMCNC+K8,2) 

1480 

C 

THE Y Yl.Xl TERMS 

1 490 


UU 160 J 3 = 1 , NC 

1500 


0=1.0 

1510 


N=CM(J3,1) 

1520 


M = CM ( J 3, 2 ) 

1530 


E=EFCP,Q,N,M) 

1540 


bETA = 6T ( N ) 

15 50 


2 3= ( 0 , 1 ,0)*N*M*PI**2/(A*b*bETAJ 

1560 

1570 


21=D6LINT(0.A,0,B,Y1) 

I f ( P . EO . N .AND. U.EU.M) THEN 

1580 


Zl=l-P*u*PI/(2*A*B))*CZl+HEM)+Z3 

1590 


ELSE 

1600 


Z1=(-P*0*PI/ (2*A*b) ) *Z1 

1610 


END IF 

1620 


Z2=DBLINT(0,A,0,b,Hl) 
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ORIGINAL PAGE Is 
OF POOR QUALITY 


page 3 


1630 

1640 

1650 

1 6b0 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
17 6 0 
1790 
1800 
1810 
1820 
1630 
1840 
1650 
1860 
1870 
1880 
1890 
1900 
1910 
1920 
1930 
1940 
1950 
1960 
1970 
1980 
1990 
2000 
2010 
2020 
203o 
2040 
2050 

2 Ob 0 
2070 
2080 
2090 
2100 
2110 
2120 
2130 
2140 
2150 
2160 
2170 
2180 
2190 
2200 
2210 
2220 
2230 
2240 
2250 
2260 
2270 
2280 
2290 
2300 
2310 
2320 
2 330 
2 340 
2 350 
2 3b 0 
2 370 
2380 
2390 
2400 
2410 
2420 
2430 


160 

C 


170 

180 


190 

200 

C 


210 

220 


230 

C 


c 


c 


c 


240 

c 


Z2=Z2*£*0/ t A * B » 2 3 
A81NC+K3, J33=Z1+Z2 
CONTINUE 

THE V VI, 11 TERNS 

DO 170 K9=1,ND 
N=CM (NC+K9 , 1 ) 

M=CM (NC + K9 , 2 ) 

D = -l .0 

E = EF (P,Q»N,M3 
bEIA=6T(N3 

23=10,-1 . 0)*(K*»2-(M*PI/B3**23/BETA 
21=DBLINTl0,A,0,b, VI 3 
IE (M.EU.O .ANO. P.LQ.N) THEN 

Zl=(K**2-(u*Pl/8)**23*lZl»R£M)/(2*Pi3*Z3 

ELSE 

Z1=IK»»2-(0*PI/63**23*Z1/C2*PI3 
END IF 

Z2=UBLINT(0,A,0,B,H23 
Z2=Z2«E»U/ 12»B*B3 
AM(NC*K8,NC+n9)=Z1+Z2 
CONTINUE 
CONTINUE 
ELSE 

DO 200 K 5 = 1 , N T 

00 190 !« 6 = 1 , NT 

AM(K5+(M-13*NT,K6+(Kl-l)*NT3=AM(K5.Nb) 

CONTINUE 
CONTINUE 
ENO IF 
ELSE 

FIND NEXT OX.OV ANO IF ALREADY SOLVED FOR LGAl) In aNO CONTINUE 
DX = DXP(K23-GXPlM 1 
UV=DYP(N23-GVP(M3 
DPT=0 

00 230 J 19 = 1 , H T P 

IF ( AbSlOX-RTC J19, 1 3 3 .LT. .OOuOl .AND. 

1 ABS(DV-RT( J19, 21 ) .LT. .00001) THEN 
DO 220 J59=1,NT 
00 210 JbO = 1 , NT 

AM((K1-1)»NT+J59, ( K 2- 13 * nT + JbO ) = AM ( (hT(j 19, 33-13 
1 »NT*J59, (RTlOl9,4)-l)*NT+J60) 

AH((K2-13»NT + J59,Ut-l3«NT + JbU) = AMl(M-l)»NT + J59, 
1 IN2-1 J *NT+obO) 

CONTINUE 
CONTINUE 
OPT= 1 . 0 
END IF 
CONTINUE 

IF NE« OX.OV SOLVE FOR COUPLING VALUES 
IF tOPT.EU.O) THEN 

IF (KIP. LT. 180) THEN 
RTP=PTP+ 1 
*T C PTP , 1 3 = (OX) 

PT l RIP , 2 1 = C DV J 
RT(RIP,3) = M 
KT(hlP,4)=K2 
ENO IF 
Trit hO» ft«5 

00 26U J b= 1 ( NC 
P=CM ( J b , 1 3 
U=CM(Jb, 23 
THE COLUMN iRS 

00 240 Jle = l ,NC 
N = CM(Jlb, 1 J 
M=CM(Jlb,2) 

V XI , X2 

E=EF(P,U,N,M) 

Z1=DBLINT(-A, A,-b,b, II 3*(K»*2-(P*PI/A3 »»23 
Z2=0bLINT(O,A,-b,B,I33 

AMUK1-1 3*NT + J6, CR2-1 J*NT*J163 = (Z1*Z23 *E 
AMHR2-1 3*NT+J6, CK1-1 3*NT+J163 = (E1*e2) *E 
CONTINUE 

V Al , V2 

DO 250 J2b=l,ND 
N = CM(J2b + NC, 1) 

M=CM ( J264NC, 2 3 
E=EF (P,0,N,M) 

Zl = CDbLlNr(-A,A,-b,b,I23 )»C-P*Q*PJ**2)/(A»8) 
Z2=DBLINT(-A,A,0,B,I4) 

AMUR1-1 3*NT+Jb, CK2-1 3 *NT+ J264NC 3 = l Z 1 ♦ Z2 3 * E 
AM UK2-1 3*NT»J6, Ul-13 *NT+ J26+NC )=( Z1 +Z2 3 *E 
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2440 

250 

CONTINUE 

2450 

2460 

2470 

260 

CONTINUE 
DO 290 L2sl.N0 
P=CM(L2»NC, 1 ) 

2480 


0 = CM ( L2*NC , 2 ) 

2490 

2500 

2510 

2520 

2530 

2540 

2550 

C 

Y Y1,X2 

00 270 L3=1.NC 
N=CM ( L3 , 1 5 
m=CM(L3,2) 

E = EF (P.U.N.M) 

Z1=(DBLINT(-A,A,-B,B,I1) ) * ( -P«Q*PI « » 2 > / ( A * 8 ) 
Z2=DBLINT(0,A,-B,B,I5) 

2560 


AM((M-lJ*NTtL2 + NC,(K2-l)*NT + L3) = (Zl + Z2l«E 

2570 


AM( (K2-1 )*NT*L2+NC, CK1-1 ) *NT*L3 ) = ( Z1 *Z2 ) *E 

2580 

2 70 

CONTINUE 

2590 

2600 

2610 

2620 

2630 

2640 

2650 

2660 

C 

Y Y 1 . Y 2 

DO 280 L4=1,ND 
H=CM(NC*L4,1) 

M=CM(nC+L4,2) 

E = EF (P.O.n.m) 

Zl=(DBLINT(-A,A.-B,b,I2))*(K**2-(0*Pl/B)**2j 

Z2=DBLlNT(-A,A,0,B,Ie) 

AM((N.l-l)*NT + L2*NC,(K2-n*NT + NC + L4)=E*(Zl + Z21 

2670 


AM((K2-l)*NT*L2 + NC,(Kl-iJ*NT + NC*L4.|=E*(Zl + Z2) 

2680 

280 

CONTINUE 

2690 

290 

CONTINUE 

2700 


END IF 

2710 


END If 

2720 

300 

CONTINUE 

2730 

310 

CONTINUE 

2740 

C 

MATRIX IS FULL 

2750 


NZ=NT*NA 

2760 

C 

NZ = TOTAL * MODES IN AN HA Y 

2770 

2760 


T Y PE* , ' PH I NT MATRIX? CY/N)' 
«CCEFT20 ,TE5T2 
IF (TEST2.E0. ' Y ' ) THEN 

2790 


2800 


DU 320 J26=1,NZ 

2810 


DU 320 J27=1,NZ 

2620 


TYPE*. 'MATklX' .J2b.J27.AMCJ2b.J27) 

2830 

320 

CONTINUE 

2840 


END IF 

2850 

28b0 

C 

FIND THE INVERSE OF THE MATRIX AND CLH V 6 VI 
CALL LEOT 1C( AM, NZ,3bO,NV,NZ, 360,0, «A, IER) 

DU 340 J 25= 1 , NZ 
V(J2S)=(0.0 J 
VI (J251=(0,0) 

2870 

2880 

2890 

330 

2900 

340 

CONTINUE 

2910 

C 

Input THE » OF DRIVEN ARRAYS, If EUUALS TOTAL » Of APERTURES 

2920 

2930 

C 

THEN ALLUM FCH EASE OF ENTRY OF PHASE SHIFTS 
TYPE*, 'INPUT THE * OF DRIVEN GUIDES' 

ACCEPT * , NDM 

2940 


2950 


10=(2.*((P1/A)**2-K**2)/E0) 

2960 


D3 = 0 

2970 


IF (NDM.E0.NA1 THEN 

2980 


IYPE*, 'CONSTANT PHASE FOR EACH? U/NJ' 

2990 

3000 


ACCEPT2U .TEST2 
IF (TEST2 . EO. ’ Y ' ) THEN 

3010 


DO 3 50 L 7 = 1 , N A 

3020 


VIC (L7-1)*NT*1 ) = ID 

3030 

350 

CONTINUE 

3040 

3050 

3060 

3070 

3080 

3090 


03 = 1 . 

ELSE 

lr (NAPE. HE. 01 THEN 

TYPE*, 'INPUT THE PhASt SHIFT bET«EEN ELEMENTS IN A R0»' 
ACCEFT * , PHG 

TYPE*, 'INPUT THE PHASE SHIFT BETWEEN RO*S' 

3100 


ACCEPT*, PhK 

3110 


NPh = 0 

3120 


03 = 1 . 

3130 
3 140 


DO 380 L71=1.NR 

PR=(L7 1-1 )*PHR*PI/180 

3150 

3160 


IF (M0DCL71 ,2) .EU.O) THEN 
DO 360 L7 2 = 1 ,NAPE 

3170 


UOO = NPH+i*( L72-1 ) * NT 

3180 


Vi (LI 00 J=ID*CEXP( (0, 1 . )*(PR+(PHG*PI/180 J *(L;2-U J ) 

3190 

360 

CONTINUE 

3200 


NPH=NPH+NAPE*NT 

3210 


ELSE 

3220 


DO 370 L73=l,NAPO 

3230 


L110=NPH*1+(L73-1)*NT 

3240 


VI (LI 10)=ID*CEXP( (0, 1 . )*(PR+(PHG*P1/ 180>*(L73-1 ) ) ) 
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3250 
3260 
3270 
3280 
3290 
3300 
3310 
3320 
3330 
3340 
3350 
3360 
3370 
3380 
3390 
3400 
3410 
3420 
3430 
3440 
3450 
3400 
3470 
3480 
3490 
3500 
3510 
3520 
3530 
3540 
3550 
3560 
3570 
3580 
3590 
3600 
3610 
3620 
3630 
3640 
3650 
3660 
367o 
3680 
3690 
3700 
3710 
3720 
3730 
3740 
3750 
3760 
3770 
3780 
3790 
3800 
3 610 
3620 
38 30 
3640 
3650 
3860 
3670 
3880 
3 690 
JQO0 
3910 
3920 
3930 
3940 
3950 
3960 
3970 
3980 
3990 
4000 
4010 
4020 
4030 
4040 
4050 


370 

380 

C 

390 

400 


410 


C 

c 


420 

43o 


440 


C 

C 


C 


C 

450 
4o 0 


470 

480 


m page is 
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CONTINUE 
NPH=NP H+NAPQ*NT 
END If 
CONTINUE 
till) If 
ENU IF 
ENO If 

If AGO NOT DRIVEN THEN INPUT I AND PHASE OF EACH DRIVEN GUIDE 
If (03. £0.0) THEN 

TYPE*, 'INPUT THE GUIDE »8S THAT ARE DRIVEN & I TPS PHASE (DEG ) ' 
DO 390 Lb=l,NDM 

ACCEPT*, VP (Lb) .PP( Lb) 

CONTINUE 

T YPE* , ' R ILL ASSUME ltlO MODE IN EACH' 

DO 400 L 8 = 1 , N DM 

V 1 ( ( VP(L8)-1 )*NT + 1 ) = 1D*CEXP( (0,1.0)*(PP(L8)*PI/180.)) 
CONTINUE 
END If 

TYPE*, 'PRINT COLUMN MATRIX? (Y/N)' 

ACCEPT20 , TEST 
If (TEST.EU. 'Y ' ) THEN 
00 410 J50=1,NT*HA 
TYPE*, ' VI : ' , VI ( J50) 

CONTINUE 
NZ=NT*NA 
END If 


USING THE RESULTS AND THE INVERTED MATRIX 
SOLVE FOR V 
DO 430 L20=1,NZ 
DO 420 L21=1.NZ 

V(l2o)=V(L20)+NV(L20,L21)*VI(L21) 

CONTINUE 

CONTINUE 

NORM 1 = ( 0 , l.)*E0/SGRT(E**2-(PI/A)**2) 

I'YpE*. 'TnE normalized output VECTOR is:' 
uU 440 LI 5 = 1 , NZ 
V15=V(L15)*N0RM1 

1 YPE*. ' V( • , LIS, ' ) = ' , V15, 1 . ' ,SGRT(kEAL( VI 5) **2*AIMAG( VI 5) **2) 
1 ATAN(A1MAG(V15)/REAL(V1S)5*180/PI 
CONTINUE 
DS = 0 


M A 1 =U 
MX2 = 0 
MX 3 = 0 
mX4 = 0 

CALCULATE THE VALLES Of ThE PATTERNS, U n 90 UEG. 
THIS DUESN'T MURK FUR A 01 MODE J t I I ! 1 1 l i 
DO 490 M10=l,91 

Th=Pl*((M10-n *2-90)/ 180 


E M 1 = 0 

Em 2 = 0 
i M 3 = 0 


EM 4 = 0 


DO 4o0 M 1 1 = 1 , N A 
DO 450 M12=1 ,NC 

£M1=EM1-JMX(TH,PI/2,CM(M12,1),CM(m12,2))*V((m 11-I)*NT 
1 ♦M12)*CtAP( (0, 1 .0)*K*SIN(TH)*DXP(MU ) ) 

EM2=EM2-JMX(Th,0,CM(Ml2. 1 ) ,CM(»12, 2) J*V( (Mll-1 ) *NT + M12 ) 
1 *CEXP((0,l.0)*K*SIN(TH)*DXP(Mll) )*COS(TH) 

E/CO 45 DEG. 

EM3=EM3-JMX(TH,PI/4,CM(M12,1),CM(M12,2))*V((m11-1)*nT 

1 ♦M12)*CCEXP((0,1.)*K*SIN(TH)*(DYP(M11)+OXP(M11))/ 

2 S0KT(2 .)))*( 1+CGS(TH ) )/2. 

Em 4=EM 4+ JMX ( Th . PI/4, CM(M12. 1),CM(M12,2))*V((M11-1)*NT 

1 +M12)*(CEXP((0,1.)*K*SIN(TH)*(DYP(M11)+DXP(M11))/ 

2 SGRT(2. ) ) ) *(COS(Th)-l ) /2 . 

THIS * AS THE E/CROSS 45 DEG. 

CONTINUE 
CONTINUE 
DO 480 M21=1,NA 

DO 470 M22=NC+1,NT 

EM3 = £N3fJMY(TH.PI/4,CMCH22,l ) , CM ( M 22 , 2 ) ) * V ( ( m2 1 - 1 ) *NT 

1 +M22)*(CEXP( (0 , 1 . ) *K*SIN(TH)*(DYP(M21)+DXP(M21) )/ 

2 3URT(2 .)))*( i-COS(TH) )/2 . 

EM4=EM4t JMY (TH , P I / 4 , CM ( M 22 , 1 ),CM(M22,2))»V((M21-l)*NT 

1 *M22)*(CEXP( (0 ,1 , ) *K*SIN(TH)*(DYP(M21 )+DXP(M21 ) )/ 

2 SGRT(2.)))*(l+C05CTH))/2. 

CONTINUE 

CuNTInUE 

Y2(M10)=SC1PT((REAL(EMI))**24(AIMAG(EMI))*«2) 
Y3(MI0)=SUKTUREAL(£M2))**2f ( A IM AG ( EM2 ) ) **2 ) 
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4060 
4070 
4080 
4090 
4100 
4110 
4120 
4130 
4140 
4150 
4160 
4170 
4180 
4190 
4200 
4210 
4220 
423o 
4240 
4250 
4260 
4270 
4280 
4290 
4300 
4310 
4320 
4330 
4340 
4350 
4 JoO 
4370 
4380 
4390 
4400 
4410 
4420 
4430 
4440 
4450 
44b0 
4470 
4480 
4490 
4500 
4510 
4520 
4530 
4540 
4550 
4560 
4570 
4580 
4590 
4600 
4b 1 0 
4620 
4630 
4o40 
4650 
4bb0 
4 b 7 0 
4680 
4690 
4700 
4710 
4720 
4730 

4740 

4750 

4760 

4770 

4780 


47 90 
4800 
4810 
4820 
4830 
4840 
4850 
4860 


490 


500 

CCC 

C 

C 


C 

C 


Y4(M10)=SORT( (REAL (EM 3 3 3 »*2* ( A IM AG (EM 3 3 3**2 3 
Y5(M10)=SGRT( (REAHEM4 3 3 **2 ♦ ( A IM AG ( £M4 3 3**2) 

AlM10) = CMo-l *2-90 
IF (MX1.LT.Y2CM10) 3 ma1=Y2(M10) 

IF (MX2.LT. Y3CW10) ) M X2 = Y 3 ( M 10 3 
IF(HX3.LT.Y4(M10)3 Mx3=Y4(M10) 

IF (MX4.LT. 15(4103) MX4=Y5(MI0) 

CONTINUE 

TYPE*. 'PLOT? ( Y/N) ' 

ACCEPT20 , ANS 
MXT=MX2 

If (MX1.GT.MX2) MXT=MX1 
DU 500 M 15=1, 91 

Y2(M15)=ObCGN( Y2(M15) ,MXT) 

Y3(M15)=0bLLN(Y3(M15),MXT) 

Y4(M153=DBCGn(Y 4(M15)»MXT) 
i5(M15)=0bCGN (Y5(M15),MXT) 

CONTINUE 

CALL AG2(X, Y2.91 ,0} 

PLOT THE PATEPNS: 

IF (ANS.EQ. ' Y 1 ) TFiEN 
PLOT TnE 4 PATTERNS! 

CALL PLOT ( X , Y 2 3 
CALL PLOT ( X , Y 3 ) 

CALL PLCJ1 ( X , Y 4 3 
CALL P LOT ( X , Y5) 

END IF 

TYPE*, 1 AN0THE8 SET OF INPUT VALUES? W/f.)' 

ACCEF120,ANS 

Al.UThER SET UF INPUT DATA f Oft THE SAME INVeRTEU MATRIX? 

IF (ANS. EG. • Y • 3 GO TO 330 

EnO 

Trit SubKOOI IN ES 
COMPLEX FUNCTIuN BT ( L 3 
REAL A, 0,>‘,L,D,tiP,ti,L 
INTEGEH P,0,N,M 

COMMON P,U,N,M,A,b,K,E,L>,OX,LY 
PI=J. 141592b 

BP=K**2-((N*PI/A)**2+(N*PI/b)**2) 

IF (bP.LT.O) THEN 

bT=(0,-l .03*SGRT(AfcS(BP) ) 

ELSE 

bT=SGRT (BP) 

ENO If 
RETURN 
ENO 

HEAL FUNCTION TP C P , N . G , P M , T 1 3 

PEAL G,HM,SP,SM,PI,G2,AL,LM,T,T1 

InTEGEh P,N 

A L= ( .001 »G) **2 

Hl=3. 141592b 

SP=F+N 

SM=P-N 

G2=2 .0*G 

0M = 0 

T=AoS(Tl } 

IF (SM.EU.03 LM= 1.0 
IF (SF.E0.03 THEN 

rp=2*(G-T)*(l-PM)/L2 

eLSe 

TP=(PM*C0S(SM*PI*T/G2) »( 1 + C0S(SP*PI 3 3 * 

1 aIN(Sf*f I*T/G2)/(SP*PI/G2)- 

2 COS(SP*PI*T/GX 3 * ( 1 + C0S(SM*PI 3 3 * 

3 (SlMSM«PI*T/G2)/(SM»HI/G2tAL3+DM*(T-G33)/G2 
END IF 

k E T U K fi 
ENO 

REAL FUNCTION k P ( P , N , G , PM , T 3 
HEAL G,PM,T,PM,P1,AL 
INTEGER P,N 
PI=J. 1415926 
A L= ( • 001 *G) **2 

9p=(G/2 . 9 3*(F'M*(SINfP*PT*T/G3-SIN(N»pT*T/f;3 ) * 

1 (l-COot(P+NJ«PI33x((P+N3*PItAL3- 

2 (SlN(P*PI*T/G)+SIN(N«Pl*T/G) J* 

3 (1-C0S((P-N)*PI))/((P-N)*PI+AL))/G 
RETURN 

ENO 

CUMPLEX FUNCTION DbL I N T ( AX , A Y , BX , B Y , F 3 
SEAL AX,AY,bX,BY,H2,S(213,X,Y,HI 

integer n 
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48 70 
4 8 8 U 
4890 
4900 
4910 
4920 
4930 
4940 
4950 
4960 
4970 
4980 
4990 
5000 
5010 
5020 
5030 
5040 
5050 
5060 
5070 
5080 
5090 
5100 
5110 
5120 
5130 
5140 
5150 
5160 
5170 
5180 
5190 
5200 
5210 
5220 
5230 
5240 
5250 
5260 
5270 
5280 
5290 
5300 
5310 
5320 
5330 
5340 
5350 
5360 
537o 
5380 
5390 
5400 
5410 
5420 
5430 
5440 
5450 
5460 
5470 
5480 
5490 
5500 
5510 
5520 
5530 
5540 
5550 
5560 
5570 
5580 
5590 
5600 
5610 
5620 
5630 
5640 
5b50 
5660 
5670 


10 

20 


COMPLEX 'F,0P 

DATA 5/1, 4, 2,4, 2, 4, 2, 4,2, 4, 2, 4, 2, 4, 2, 4, 2, 4,2, 4,1/ 
h = 20 

H1=(AY-AX)/N 
h2=Idi-6X)/N 
DP = 0 

CO 20 11=0, N 
DO 10 1.2 = 0 , N 
A=Ll»hlfAX 
Y=L2»H2+6X 

DP=DP+FCX, Y1*S(L1+1 J*S(L2+1) 

CONTINUE 

CONTINUE 

D6LIM = OP»H1*H2/9.0 

RETURN 

£N0 

COMPLEX FUNCTION THU, 8) 

SEAL U , « , A , b , E , K , 0 , AL , TP 

COMPLEX lP , 

integer p,u,n,m 

CGMMUN P.U,N,M,A.6,K,E,D,DX,DY 
AL=( A*.001)»»2 

UsCeaPC (0, -1.0) *K*SUPT( U* *2 + 8 4*2) 1 /SORT C U* *2 + »* * 2 ♦ AL ) 
ip=ip»tp(P,n,a,d,u)8TPco,«,b,-;.o*d,h)8E 
IF (P.EQ.N .AND. Q.EO.MJ Then 
11=1P-4.0/SOHTCU*«24«**2+AL) 

ELSE 
Y 1 = YP 
ENO If 
Rtl'URN 
End 

Cu-PLEX FUNCTION H 1 C X , <• 3 
REAL 0, A , 6 ,K , E, X , » , AL, TP, PI 
CUMpLt-A HP 
INTEGER P,G,N,M 

COMMON P.Q,N,m,A,B,K,E,D,LX,UY 
AL=(A».O0l ) »»2 
PI=3. 1415926 

HP=CEXP((0,-1.0)«N«SGRT(8**2tX*»2))/SORT(X»»2+»**2fAL) 

1 -CCS(P»PI) *CEXP( (0,-1.0)»K*SCihTllA-X)**2*<i**2))/ 

2 SahHCA-aJ** 2 + 8**24AL) 
nP=nP«SINCN*FI»X/A) 
iU=HP«TP(G,M,fa,-1.0,«) 

H El URN 
END 

COMPLEX FUNCTION H2(U,Y) 

MEAL U,Y,L,K,A,B,E,AL,TP,Pl 
IhTEllK p,g,n,m 
COMPLEX np 

COMMON P,G,N,M,A,B,N,E,D,DX,0Y 
AL=( A*.u01 ) * » 2 
Pl=3. 1415926 

hP=CEXP( 10 ,-l , 0 1 *R*SS2PT(U«*2 + Y**2) )/SGRT(U»*2+Y**2«-ALH 

1 C0S(0*PI) «CEXP( 10,-1 .01 *a*S0PTCU»*2+(6-Y)*«2 11/ 

2 SuRT(0«» 2* (B-Y 1 »=2+ALl 
HP=hP*SlN(M*PI*Y/6) 
h2=HP*TP(P,N,A,-1.0,U) 

R t TO P n 
tNU 

COMPLEX FUNCTION G(X,Y1 

HEAL X,Y,P1,K,R,AL 

COMMON P,0,N,M,A,8,K,E,D,DX,DY 

PI=3. 1415926 

R=SoRT(X«*2+Y**2) 

G=C EXP U 0,-1 .01*n*R)/(2*PI*K) 

REIuHn 

ENO 

COMPLEX FUNCTION I|(U,»l 
HEAL U,H,TP,RP,A,b,K,D,LX,DY 

complex g,ip 

INTEGER P,o,N,M 

COMMON P.Q,N,M,A,B,K,E,D,DX,LY 

:? = f 1/4. 5) *(T? P ,N , A , I . J ,L i t R P 1 P , N , A , 1 .3, J) ) * 

1 CTP(G,M.b,-1.0,8j+HPCO,M,b,-1.0,»)l« 

2 G( U-OX , »-DY ) 

1 1 = I P 

RETURN 

iO 


COMPLEX FUNCTION I21U,«) 

REAL U,«,I?,rtP,A,d,A,j,OX,DY 
COMPLEX G, IP 
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5680 

5690 

5700 

5710 

5720 

5730 

5740 

5750 

5760 

5770 

5760 

5790 

5600 

5810 

5620 

5630 

5640 

5850 

5860 

5870 

5660 

5890 

5900 

5910 

5920 


5930 
5940 
5950 
5960 
5970 
5980 
5990 
6000 
6010 
6020 
6030 
6040 
6050 
bOoO 
6070 
6080 
6090 
6100 
6110 
6120 
6130 
6140 
o 1 5 0 
6160 
6170 
6180 
6190 
6200 
6210 
6220 
6230 
6240 
6250 
6260 
6270 
6280 
6290 
6300 
6310 
6320 
6330 
6340 
6350 
6360 
6370 
6380 
6390 
6400 
6410 
6420 
6430 
6440 
6450 
6460 
6470 
6480 


INTEGER P,u,N,M 

COMMON P.G,N.M,A,8,K,E,D,DX.bY 
IP=(l/4.6)*(TP(P,h,A,-I.O,U3*SP(P#N,A,-1.0,U})* 

1 (IPCO, M, 6,1.0, • )+kPCQ,M, b, 1.0,8))* 

2 GCU-OX, »-DY ) 

12 = I P 
RETURN 
Eft 0 

COMPLEX J UNCTION I 3 C X , « J 
HEAL X,M,TP ( KP,A ( b,K ( D,OX,l/Y,PI,E 
COMPLEX G,IP 
INTEGER P,G,N,M 

COMMON P,G,h,M, A,8,K,E,D,DX,0Y 
81=3.1415926 

IF=IP*PI/(A«A*2) )*S1N(N*PI*X/A3 

IP = IP*CG(0-X-UX ,*-DY3-C0S(P*PI 3*G( A-X-DX, a-DY 3 ) 

13 = IP*CTP(0,K,B,-1 .O,m3 + KP(0,M,d,-1 .0, « 3 3 
HETURN 

f NO 

COMPLEX FUNCTION I 5 l X , w ) 

REAL A»«/1R.HP/A ( B,K,0»DX / 0Y»PI,E 

complex G,1P 

INTEGER P,Q,N,M 

CUMMON P,G . N ,M , A , 8 ,K , E,U,CX, DY 
PI=3. 1415926 

lP=(G»PI/(a»A»2))»SINCN»PI*X/A) 

IP = IP*(G(0-X-0X, V.-OY )-CUSCP»PI J*G( A-X-UX, *-DY ) J 
i5=lP4(TP(G # M,b,-1.0,«)+RPlo,M,B,-1.0,«.)) 

RETURN 

EftD 

COMPLEX FUNCTION I4(U,Y) 

HEAL U.Y , T F , p, P , A,B,K,U,DX,Q'i ,P1,E 
COMPLEX G , I P 
INTEGER f,u,n,m 

COMMON P,U,N,M,A,B ( K,E»D#DA / DY 
HI=3. 141592b 

lP=(P*PI/(2.0*A*B))*SIN(M*Pl*Y/b) 
1P=1P*<G(U-UX,0-Y-UY3-CCSCG*PI3*GCU-DX,B-Y-UY 3 3 
I4=IP*(TP(P,N,A,-1.0,UJ*RP(P,N,A,-1.U,U33 
RETURN 

end 

COMPLEX JUNCTION Ib(U,Y) 

HEAL U.Y.TP,PP,A,B,K,D,DX,DY,PI,E 
COMPLEX G,IP 
INTEGER P,U,N,M 

Common P,O.N,M,A,b,K,E,U,OX,UY 
HI=3. 1415926 

IF=(U*PI/(2.U*B*633*SINCM*PI*Y/b3 

iP=IP«(G(U-DX,0-Y-DY3-C0SC0«P13*GCU-C'X,c-Y-LY33 

Ie=lP*irP(P,N,A,-I.U,U3+KP(P,N,A,-1.0,U)) 

RETURN 

EftO 

kEAL JUNCTION efcp,o,n,m; 
real, E 

integer p , u > n , m 

E= 1 6 . 0 

IJ lG.EU.03 E=E/2.0 
IF CM.EO.O) E=E/2.0 
IF CP.E0.03 E=E/2.0 
IF CN.EG.03 E = t/2 . 0 
eF = SGR TIE) 

RETURN 

END 

COMPLEX JUNCTION JMX (T.PH.NI ,M1 3 
PEAL I.Pri.El ,PI,KX,RY,K, A,P,G1 
I n TE gER N 1 , M 1 

COMMON P,0,N,M,A,6 r K,E,D#DX,LY 
P I =3 . 141592b 
G 1 = ( ,001«A3**2 
E 1 =4 , 0 

IK Chi. £0.03 E 1 =E1 / 2 
IF ( M 1 . tU . 0 3 El = El/2 
E 1 = SURT C El 3 
KX=6iN(T3»COSCPH3»K 
nY=SINCT 3*SIN(PR3*K 
Ir (PH.EU.03 THEN 
IF (Ml. £0.03 THEN 

JMX=tl« (NI*° 1 3 «( l-CGSC Nl *PI 3 *CEXPC ’ .'J;«k.X*A) 3 / ( (Nl *01/ A 

1 **2-KX**23 
ELBE 
JMX=0 
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PAGE is 
OF POOR QUALITY 


PAGE 9 


6490 

6500 

6510 

6520 

6530 

6540 

6550 

6560 

6570 

6580 

6590 

6600 

6610 

6620 

66 30 
6640 
6650 
b6eO 
667 0 
06 8 0 
6690 
6700 

67 10 
6720 
o 7 30 
6740 
e750 
6 7b0 
6770 
6760 
6790 
6800 
6810 
6820 
6830 
6840 
6850 
6860 
6870 
6660 
6890 
6900 
6910 
6920 
6930 
6940 
6950 
b9e0 
6970 
6960 
e990 
7000 
7010 
7020 
7030 
7040 


C 

C 


C 


ENG IF 
ELSE 

IE (M1.EG.0 .AND. T.ECl.O) THEN 

JHAS-El * (1-C0S(N1 *PI 3 3 *B/ (N1*PI/ A) 

ELSE 

JMA=E1*(N1*PI/A) *( (0, 1 ,0)*M)*( 1 -COS(N 1»P1) * 

1 CEAPC (0. 1.0)*KA*A3 )*<CEAPl CO, 1 . 0 ) »A Y • B ) *C0S ( m 1 » P 1) - 1 ) / 

2 ( UN1»PI/A)**2-KX**2+G1 J»( C“1*PI/B)*«2-K***24G1J ) 

END IF 

END IF 
RETURN 
END 

ADD SYMMETRIC MATRIX PART 

A i/O EAoIER PHASE ENTRY 

COMPLEX FUNCTION JMY CT.PH.Nl ,M1 ) 

REAL T,Ppi,E1.RI,KX,nY,K,A,B,G 1 
INTEGER N 1 .Ml 

COMMON P,G,N,M,A,8,K,t,D,DX,DY 
PI=3. 1415920 
Gl= t .001 »A) *«2 
El=4. 

IKH1.EU.0) El=El/2. 

II (Ml.Eo.01 E 1 =E 1 / 2 . 

E 1 =SOKT (Ell 
rX=K»SIn(T) «COS<PH) 
rY=6»S1N(T)*S1N(Pm) 

IF (T.EO.O) ThEN 
jMY = 0 
ELSE 

JMYSE1«(M1»PI/B)«( 1-CLS(M«PI )*CEXP( (0, 1 . )*M»A) )* 

1 lO , 1 . ) *KX* (C0S(M1«PI )*CEXPUO, 1 . ) *AY«6) - n/ ( ( 

2 IN1«PI/A)**2-rX*» 2)*( (M1»PI/B)»*2-NY**2)+GU 
END 1 r 

RETURN 

EkD 

SUBROUTINE PL('TCX.Y) 

KFAL X C 9 1 J , Y (91) 

CAuL I N 1 XT ( 240 1 
CALo oINITT 
CAUL UblMY 1-60. ,0J 
CAuL h P T S ( 9 1 } 

CALL CHECKX.Y) 

CmLL DSPLAYU, f) 

CAUL IKPAoSE 
CalL IN' ITTC 240 ) 

CAUL f 1MTTC0 , 767 ) 

RETURN 

END 

This 1NCLUDFS THE 45 UE&. PLOTS 
READ JUNCTION U8CUNCX.Y) 
reAL X , Y 
X = A/Y 

IF CX.LE. 0.001) Xs.001 
ubCUis = 2o»D0G10(X) 

RETORN 

ENu 


APPENDIX B 


USER INSTRUCTIONS 

The computer program discussed in this report and presented in 
Appendix A is currently on the Research Computer Labs VAX under the 
file name: MCOUPL. It is written in FORTRAN and can handle large 
matrices, on the order of 300 x 300 (the user is warned that tests of 
matrices larger than 125 x 125 have never been done and round-off 
error may be a problem). The following is the user instructions, and 
before proceeding the notation used must be explained. In the exam- 
ples shown, the underlined terms are the user entered responses, while 
the rest are the computer generated phrases and questions. The * 
signifies a return. The examples to be discussed are presented at 
the end of this section, with added circled numbers signifying the 
step numbers referred to in the text. 

The program is written so that the computer will prompt the user 
for the necessary data. There are a few factors that are held constant 
throughout the array and they are entered first. The constant factors 
are the a and b dimensions of the apertures entered in wavelengths, 
the spacing between each row and between each element in each row again 
in wavelengths, the wavelength entered in meters and finally the num- 
ber of and type of ib and $ to be used in the approximation. 

nm nra 

These are values that typically remain constant throughout standard 
waveguide arrays. In the sample runs, shown at the end of this section, 
steps 1-3 are the ones where the above mentioned data is entered. In 
step 1, after inputting the number of 4 > and $ the computer will ask 
for the mode index nm one at a time. The order that these are entered 
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must always be the first modal function entered. Steps 2 and 3 are 
self explanatory. 

Step 4 is the point where the user has the choice of whether to 
take advantage of an array configuration entry system that is designed 
specifically for arrays of rectangular or triangular lattice or to use 
the general entry system. In the steps that will follow, the user 
enters the array configuration and if it is one of the regular lattices 
it can be entered much easier under the system mentioned above. As we 
will see later, this approach will also allow for easier entry of 
regular phase shifts in the driven fields used when scanning the beam. 
The regular lattice approach will be looked at first and the general 
method will be discussed in a later example (see steps 12 and 13) 

The array must have at least 2 rows to use the regular lattice method. 

After a "Y" is entered (note: in all the Y/N questions, the Y must 
be upper case) the computer is in the regular lattice mode. Under this 
system the user enters the number of apertures and the row off- 
set (in X) for the even numbered rows and for the odd numbered rows 
once each. The way the computer numbers rows and apertures is shown 
in Fig. A. 1 The rows always count up and the apertures from left to 
right. The configuration entered in step 4 is the one shown in Fig. 

A. la. Figure A. lb shows what is meant by a row off-set and as to how 
it allows for entry of triangular lattice arrays. If the array of 
Fig. A. lb were entered instead, 3 sows would have been entered in 
step 3 and for the odd numbered rows an offset of .5 would have been 
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handled by each configuration entry system. 

After this, the computer will print out the D and D of each 

x y 

aperture in meters in the order of the aperture numbers. The first 

number is the aperture number, the next the D and the last the D 

x y • 

After this the computer will calculate the admittance parameters. 

Once this is complete, the user has the opportunity to view them by 
entering a "Y" in step 5. The computer will then print first the 
column number and then the row number of each parameter's location in 
the admittance matrix. The value of the admittance is given in 
standard FORTRAN complex number representation (the real part comes 
first, followed by the imaginary part, seperated by a comma). If 
a "N" is entered the computer would go directly to Step 6. 

Step 6 is where the user tells the computer how to drive the 
array in calculating the modal amplitudes and then the radiation 
patterns. The user must enter how many elements, and if it's less 
than the total number which elements are to be driven. Also the phase 
of the driving field of the elements must be entered. The computer 
alwavs drives the first i b of each driven aperture with a TE n 
mode, so the user must make sure to always enter the 10 mode first in 
step 1 as was done here. Under the regular lattice system the user has 
3 different ways to enter the driven element data. The first is shown 
in steps 6 and 7. If the number of driven elements equals the total 
number of elements then the computer will ask if a constant phase for 
each should be used. If "Y" is entered here, then all will be 
equally driven with zero phase. The case of "N" entered under con- 
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stant phase is shown in step 11. Here the computer asks for the phase 
shift between each element in a row (call it g here) and the phase 
shift between each row (call it a here), both entered in degrees. 

The added phase for each aperture then equals (n-l)g + (m-l)a where 
n is the column number and m is the row number of the aperture. This 
way, aperture #1 will always have 0° of added phase shift. For 
rectangular lattices this allows the user to easily scan the beam. 

There is a problem with this system when using a triangular lattice. 

Since the computer always assigns the first aperture of each row to 
column irl and since the apertures aren't co-linear in the vertical 
direction this will cause problems in these cases. This will be 
corrected in the next update of the software, but for now if you wish 
to scan beams in a triangular lattice use the general approach and 
enter the phase of each aperture one at a time. The third and final 
situation is shown in step 10. If less than the total number of 
elements is entered then the user must enter the element number and it's 
added phase for each element to be driven. After each element number 
and it's corresponding phase, a return must be entered as shown 1 ere or 
else the computer will not input the data properly. The steps tha ' follow 
each of these situations are the same and will be discussed next, using 
the example steps 7-10. 

Shown after step 7 is the next question: That of whether to view 

the tested incident magnetic field column matrix (see page 27.) If 
a "Y" is entered, the elements are printed out in order, as shown. 
Notice how it's the first modal function of each aperture that is 
driven and with zero phase in this example. Directly after this, or 
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after a "N" is entered above, the computer prints out the normalized 
modal amplitudes of each modal function. The order of these amplitudes 
is the order of the apertures (as numbered before), where within each 
aperture the order of the modal functions is the same order as entered 
in Step 1. In this example, there were 4 modes per aperture so the 
first 4 amplitudes are the ones of aperture //I , the next 4 for aperture 
If 2 and so on. The modes of each aperture then are in the order of 

^10’ ^ 12 ' ^ 30 ’ anc ^ ^i? as enCere d in Step If 1. So the first ampli- 
tude is the •(> of aperture If 1 , the second 1 i> and so on. 

In Step 9, the user decides whether to view the plots of the 
radiation pattern. If a "Y" is entered, then the patterns will 
appear one at a time in the following order. 1) $ = 90° co-polarized , 

2) (j> = 0° co-polarized, 3) <j> = 45° co-polarized and, 4) <j> = 45° cross- 

polarized. After each plot has been viewed, hitting the return will 
bring up the next pattern until the last plot. Hitting the return 
after viewing the last plot brings back the program. 

After the plot question, the user has a choice, to input a new 
set of values for the driven elements or to exit the program (see 
Step 10). A "N" entered here and the program exists. A "Y" and the 
computer returns to Step 6 and proceeds exactly as before. In the example 
shown, a "Y" was entered in Steps 10, and right before Step 11 in 
order to enter different sets of driven elements for the same array. 

This way the user can scan the beam or change the number of driven 
elements several times on the same array configuration without going 
through the redundant procedure of having the computer recompute the 
admittance parameters. 
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Finally Steps 12 and 13 show how the general array configuration 
entry system works. If a "N" was entered when the computer asked 
if the regular lattice approach was to be used, the user will be promp- 
ted to input the number of elements and row offset for each row in the 
array one at a time. The example shown in 12 and 13 is the array of 
Figure A. lb. After these steps the computer proceeds exactly as before 
with one exception. When the user drives all the elements as in Step 
6 he has the choice to drive all with constant phase, as before, but 
if a "N" is entered there the user then must enter the element num- 
bers and their phases one at a time similar to what was done in Step 
10. The user does not have the option to use the phase method shown 
in Step 11. If less then the total number are driven, it is handled 
exactly as before. 

As a final comment, the meaning of the normalized modal ampli- 
tudes will be discussed. The printed values are the normalized equi- 
valent magnetic current amplitudes (J^) , not the normalized aperture 
field amplitudes (see Appendix C for the differences between the two). 
The currents are always normalized by an incident term of 0° 

added phase, so when entering the phases of the driven modes always 
make the reference phase 0°. The usefulness of this type of result 
is readily apparent (see App. C). If only aperture i is driven (with 
0° phase), then the resulting values of the normalized i|j^q amplitude 

term of aperture j is the S . term. For arrays of 

10 1 10 J 

square apertures the value of the normalized medal amplitude of the 

0 n term in aperture j then is the S „ . , T _ . • Finally, the nor- 

Ui it 01 3 

malized value of the driven aperture's amplitude gives 1 4- I\ for 
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that guide. So by driving only one guide with zero phase the computer 
will output the scattering parameters for coupling between the driven 
guide and it's neighboring undriven guides and the reflection co- 
efficient of the driven aperture. 
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Figuie A.l.i-d Different Arrays and Which Kntry System Can Be Used. 
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INPUT- THE ♦ PSI AND PHI FUNCTIONS 

AJL* 

> INPUT THE 1PSI COEFF.: 

AjJL* 

INPUT THE 2PSI COEFF.: 

1,2 * 

INPUT THE 3PSI COEFF.: 

3.0 * 

INPUT THE 1PHI COEFF.: 

1 * 2 .* ■ 

INPUT A /LAMBDA .B/LAMBDA AND LAMBDA 

. 6 . . 6 . . 01 5 * 

ALL SPACINGS ABE TO BE GIVEN IN F PACTIONS OF LAMBDA 
INPUT THE # OF ROWS AND SPACING BETWEEN THE BOWS 
2 1 .* 

INPUT THE SPACING BETWEEN ELEMENTS WITHIN A ROW 
1 , * 

PECULAP LATTICE ARRAY? (Y/N) 

Y * 

INPUT THE EVEN ROWS * OF GUIDES & ROW OFFSET 


© 

© 

© 

© 


3.0 * 

INPUT THE ODD ROWS # OF GUIDES & ROW OFFSET 


-UL* 

DX ,DY 
DX ,DY 
DX ,DY 
DX.DY 
DX.DY 
DX.DY 

PRINT MATPIX? 

X* 

MATPIX 

MATPIX 

MATPIX 

MATRIX 

MATHIX 


1 0. OOOOOOOE+OO 

2 1. 5000000E-02 

3 2. 9999999E-02 

4 0. OOOOOOOE+OO 

5 1. 5000000E-02 

6 2. 9999999E-02 
(Y/N) 


0. OOOOOOOE+OO 
0. OOOOOOOE+OO 
0. OOOOOOOE+OO 
1 . S000000E-02 
1 . 5000000E-02 



1 1 (12.27046966552734,-517.9987182617186) 
1 2 C.-30. 51866722106934 ,54. 47422790527344) 
1 3 (53.73537826536086,-70.30246734619141) 
1 4 (71.57903289794922,-26.94529342651367) 
1 5 (1.308216094970703,19.72079849243164) 


INPUT THE # OF DRIVEN GUIDES 
6 

CONSTANT PHASE FOP EACH? (Y/N) 

Y + 

PRINT COLUMN MATPIX? *(Y/N) 

X.* 

VI: ( -1.2 110 38 07718 35 904 E + 16 , 0 . 0000000000000 000E + 00) 
VI : (0. OOOOOOOOOOOOOOOOE + QO ,0. OOOOOOOOOOOOOOOOE + OQ) 
VI: (0. COOOOOOOOOOOOOOOE+O0,0. 0 000 00000000 0000E+00) 
VI: (0. 000000000000 000 0E + 00 ,0.0 000 0000 OOOOOOOOE + 00) 
VI : (-1.2 11 03807 71835904 E+16, 0.000 000000 OOOOOOOE+OO) 
VI: (0. 0000 000 00000 000 OE+00, 0.0 00000000000 0000E+00) 
VI : (0. 00000 0000000 00 00E+ 00 ,0. 0000 OOOOOOOOOOOOE+OO) 
VI : (0. 00 0000000 OOOOOOOE + OO ,0. OOOOOOOOOOOOOOOOE+OQ) 
VI: C-1- 2 1103807718 35904 E+16, 0 . OOOOOOOOOOOOOOOOE+OO) 
VI : (0. COOOOOOOOOOCCOOOE + CO ,0. OOOOOOOOOOOOOOOOE+OO) 
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ORIGINAL 
Of POOR 


PAGE IS 
QUALITY 


THE NORMALIZED OUTPUT VECTOR IS: 

V( 1)= (0. 9525498,-6. 508381 7E -02) 

V( 2) = (0. 1387619 , -9. 7904637E-02) 

V( 3) = (2.6979776E-02 ,-4. 0551741E 

V( 4 ) = (-2. 4122884E-02 ,6 .4293578E 

V( 5 ) = (0.9735299 , -4. 6621375E-02) 

V( 6 ) = (0.1327942,-9. 0362079E-02) 

V( 7 ) = (2. 4720680E-02.-4.1649926E 

V< 8 ) = (-1.8278474E-02.7.0564717E 

V( 9 ) = (0.9525498 , -6. 5083809E-02) 

V( 10)= (0.1387620,-9. 7904630E-02) 

V( 11)= (2.6979776E-02,-4.0551744E 

V( 12)= (-2.4122884E-02,6.4293578E 

V( 13)= (0.9525498, -6. 5083809E-02) 

V( 14)= (0. 1387620,-9. 7904630E-02) 

V( 15)= (2.6979776E-02,-4.0551744E 

V( 16)= (-2. 4122884E-02 ,6.4293578E 

V( 17)= (0.9735299 , -4. 6621375E-02) 

V( 18)= (0.1327942 , -9. 036207SE-02) 

V( 19)= (2. 4720680E-02,-4.1649926E 

V( 20)= (-1.8278474E-02, 7. 056471 7£ 

V( 21)= (0.9525498 , -6. 5083817E-02) 

V( 22)= (0.1387619 , -9. 7904637E-02) 

V( 23)= (2.6979776E-02,-4.055174lE 

V< 24)= (-2.4122884E-02,6.4293578E 

PLOT? (Y/N) 


0.9547706 
0.1698 24 0- , - 

02), 4.8706800E 

■02), 6.867C064E 

0.9746456 
0.1606225 
■02), 4 . 84 33 751 E 

02) , 7.2893634E 

0.9547706 
0.1698240 
02), 4. 8706800E 

02), 6.86 70064E 

0.9547706 
0.1698240 
02), 4. 67C6600E 

02), 6 . 86 7006 4E 

0.9746456 
0.1606225 
02) , 4.84 33751E 

02), 7. 2693634E 

0.9547706 
0.1698240 
02) , 4. 8706600E 

02) , 6.86 70064E 


-3.908711 
-35.20520 
■02 -56.36349 

•02 -69.43357 

-2. 74)743 
-34. 23400 
-02 -59.30940 

-02 -75.47779 

-3.90671C 
-35.20520 
-02 -56.36349 

■02 -69.43397 

-3.908710 
-35.20520 
-02 -56.36349 

-02 -69.43397 

-2. 741743 
-34. 23400 
-02 -59.30940 

-02 -75.4 777 9 

-3.908711 
-35.20520 
•02 -56.36349 

-02 -69.43397 


ANOTHER SET OF 
y * 

INPUT THE * OF 


INPUT VALUES? (Y/N) 
DRIVEN GUIDES 


INPUT THE GUIDE i' S THAT ARE DRIVEN SIT'S PHASE (DEG) 

2 . 0 * 

J.Q: 

WILL ASSUME TE 10 MODE IN EACH 
PRINT COLUMN MATRIX? (Y/N) 

N « 

THE NORMALIZED OUTPUT VECTOR IS: 

(0.9562839,-2. 7126919E-02) , 


V( 

1)« 

V( 

2) = 

V( 

3) = 

V( 

4 ) = 


9566686 

1488423 


(-1. 324 4 004E-0 3, 5. 052013 3E-02) , 


-1.624875 
-35.13311 
4.48 50983E-02 -59.28555 

5. 0537493E-02 -88.49832 


ANOTHER SET OF INPUT VALUES? (Y/N) 

X * 

INPUT THE t OF DRIVEN GUIDES 
6 * 

CONSTANT PHASE FOR EACH? (Y/N) 

_N_* 

INPUT THE PHASE SHIFT BETWEEN ELEMENTS IN A ROW 

JO* 

INPUT THE PHASE SHIFT BETWEEN ROWS 
30 * 

PRINT COLUMN MATRIX? (Y/N) 



ALL SPACINGS ARE TO BE GIVEN IN FRACTIONS OF LAMBDA 
INPUT THE # OF ROWS AND SPACING BETWEEN THE ROWS 

■ 3 , 1 * 

INPUT THE SPACING BETWEEN ELEfrENTS WITHIN A ROW 


1 * 


REGULAR LATTICE 

JL* 


INPUT 

3.0* 

THE 

# 

OF 

INPUT 

3-0* 

THE 

# 

OF 

INPUT 

THE 

# 

OF 


3.0* 

DX.DY 

DX.DY 

DX.DY 

DX.DY 

DX.DY 

DX.DY 

DX.DY 

DX.DY 

DX.DY 


ARRAY? (Y/N) 

APP. IN ROW# 1 

APP. IN ROW# 2 

APP. IN ROW# 3 

1 O.OOOOOOOE+OO 

2 1. 5000000E-02 

3 2. 9999999E-02 

4 O.OOOOOOOE+OO 

5 1 . 5000000E-0 2 

6 2. 9999999E-02 

7 O.OOOOOOOE+OO 

8 1. S000000E -02 

9 2.9999999E-02 


AND THE ROW OFFSET 

AND THE ROW OFFSET 

AND THE ROW OFFSET 

O.OOOOOOOE+OO 
O.OOOOOOOE+OO 
O.OOOOOOOE+OO 
1.5000000E-02 
1. 5000000E-02 
1. 5000000E-02 
2. 9999999E-02 
2. 9999999E-02 
2.9999999E-02 


APPENDIX C 


RESULT NORMALIZATION 

The results shown in Chap. 3 are the normalized aperture electric 

field modal amplitudes, normalized by the incident field amplitude. 

The program, on the other hand, will output the normalized magnetic 

current (J ) modal amplitudes. The differences between the two and 
— m 

the reason for outputting the normalized currents is the topic of this 
appendix. 

The computer model solves for the J modal amplitudes. Since 

m 

E = a x J and J = y 1 C ^ a+V'D^a 

— apert. — z — m — m < nm nm — x 4_< nm nm — y 

n ’ m n,m 

then E =-y^Ctpa+y'D$ a. 

—apert . nm nm — y L-> nm nm — x 

n,m n,m 

INC 

The E 4 " 1 normalization term has an extra minus sign relative to the J 
— -m 

normalization term. The difference between the normalized J and the 

— m 

normalized E is that the signs on the D terms are changed and the 
— a nm 

a and a 's interchanged. Thus taking the negative of the 0 
— x — y nm 

amplitude terms from the computer outputs and assigning the proper a^and 

a^ is required to get the normalized _E^ functions. 

The reason for outputting the results in terms of normalized J 

— m 

is that this leads directly to the scattering parameters. Scattering 
matrix parameters are defined using the TE waveguide modal functions. 
These modal functions by definition have opposite signs for the _a 
and _a terms (see the modal function used by Harrington in Sec. 2.3 
for example). This negative sign cancels the one needed to convert 
the currents to field quantities. So by leaving it this way, if only 
one aperture is driven in an array of square guides, the value of the 
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normalized amplitude in a neighboring guide gives the scattering 

parameter for coupling between the TE ^ of the driven guide and the 

TE 01 of the undriven guide. If E terms were outputted it would 
01 — a 

equal the negative of this S parameter. The value of the normalized 
amplitudes are the scattering parameters between the driven TE^q 
and the undriven TE ]_q in both cases. It must be remembered that 
this is only so if only one guide is driven in the array with zero 
degrees added phase. 

The value of the amplitude used to normalize the E term is: 

6 A ~ a 

„INC P 10 H 10 
E = - — . 



